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ABSTRACT 

Using combinations of Ha, ultraviolet (UV), and infrared (IR) emission, we estimate the star for- 
mation rate (SFR) surface density, Ssfr, at 1 kpc resolution for 30 disk galaxies that are targets 
of the IRAM HERACLES CO survey. We present a new physically- motivated IR spectral energy 
distribution-based approach to account for possible contributions to 24/im emission not associated 
with recent star formation. Considering a variety of "reference" SFRs from the literature, we revisit 
the calibration of the 24/im term in hybrid (UV+IR or Ha+IR) tracers. We show that the overall 
calibration of this term remains uncertain at the factor of two level because of the lack of wide-field, 
robust reference SFR estimates. Within this uncertainty, published calibrations represent a reason- 
able starting point for 1 kpc-wide areas of star-forming disk galaxies but we re-derive and refine the 
calibration of the IR term in these tracers to match our resolution and approach to 24/im emission. 
We compare a large suite of Ssfr estimates and find that above Ssfr ^ 10~^ Mq yr~^ kpc~^ the 
systematic differences among tracers are less than a factor of two across two orders of magnitude 
dynamic range. We caution that methodology and data both become serious issues below this level. 
We note from simple model considerations that focusing on a part of a galaxy dominated by a single 
stellar population the intrinsic uncertainty in Ha and FUV-based SFRs are '^0.3 and '^0.5 dex. 
Subject headings: galaxies: evolution — galaxies: ISM — ISM: dust, extinction — stars: formation 



1. INTRODUCTION 

Comparing the distributions of recent star formation 
and interstellar gas can illuminate the preconditions for 
star formation, the impact of local conditions on the con- 
version of gas to stars, and the process of galaxy as- 
sembly. Making these comparisons quantitative requires 
local estimates of the recent star formation rate (SFR) 
that span whole galaxies. Such estimates can be re- 
constructed from the observational signatures of massive 
stars, which live only a short time after they are born. 
These signatures include recombination line and ther- 
mal radio emission, both produced by ionizing photons 
emitted by the most massive stars. In regions of active 
star formation, direct photospheric emission from mas- 
sive stars also dominates ultraviolet (UV) emission and 
provides the bulk of starlight reprocessed by dust into 
infrared (IR) emission. 

This paper investigates the use of such observations 
to estimate the surface density of recent star formation, 
^SFR, at 1 kpc resolution in nearby galaxies. We focus 
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on 30 galaxies that are targets of the IRAIV0 HERA- 
CLES Large Program ( Leroy et al. 2009). HERACLES 
surveyed CO J = 2 ^ 1 emission and we have H I maps 
for each targe t, mostly from the VLA THINGS survey 
([Walter et al.| [2008). As a result, we know the distribu- 
tions of both the molecular (H2) and atomic (H I) inter- 
stehar medium (ISM). THINGS and HERACLES aim to 
constrain theories of star formation via comparisons of 
HI, H2, SFR, and environment at matched resolution 
across large areas. Ssfr estimates matched to the gas 
maps play a key role in such tests a main goal of this 
paper is to motivate our approach to such comparisons. 

There is no definitive approach to estimate Ssfr at kpc 
resolutions. Indeed, less attention has been paid to esti- 
mating HsFR for large parts of galaxies than to SFRs for 
individual regions or whole galaxies. We therefore begin 
at a low level and focus on issues related to this interme- 
diate scale as we work towards Ssfr estimates. These 
include the scatter in SFR estimates that arises from 
considering individual stellar populations, the fraction 
of IR emission not associated with recent star formation, 
and the appropriate calibration for "hybrid" UV+IR or 
Ha+IR tracers. 

We begin by defining our approach, which follows other 
recent work by hybridizing UV or Ha emission with IR 
emission (Sections [3] and [4j). We consider how UV and 
Ha-based estimates will be affected by isolating indi- 
vidual stellar populations rather than considering an in- 
tegrated, continuously star-forming population (Section 
[3]). This will be an increasing issue moving from inte- 
grated galaxi es to sub-galactic sc ales and we examine 
Starburst99 ([Leitherer et al.l |1999[ ) simulations to esti- 
mate the magnitude of the issue. We discuss basic trends 

9 IRAM is supported by CNRS/INSU (France), the MPG (Ger- 
many) and the IGN (Spain). 
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among observed intensities of UV, Ha, and IR emission 
at 1 kpc resolution, including their correlations and rela- 
tive magnitudes (Section |4]). We then explore the effects 
and potential magnitude of IR emission not directly asso- 
ciated with star formation (Section[5]). We adopt a physi- 
cally motivated approach to this topic, interpreting 24/im 
emission through the lens of dust models fit to the whole 
IR spectral energy distribution (SED) and test its effect 
(Section [5j). Based on these calculations, we explore the 
calibration of empirical "hybrid" (IR+Ha and IR+UV) 
tracers for our regime (Section [6]) and discuss previous 
work on this topic. We emphasize the substantial un- 
certainty that still remains in the absolute calibration of 
such tracers despite much recent attention. We conclude 
by comparing a large suite of Ssfr estimates (Section 
[7j) and then discussing uncertainties in SFR estimates, 
outlining our recommended approach to estimate Hsfr, 
and noting several implications of our work for studying 
star- formation in galaxies (Section [8|). We present our 
conclusions in Section [9l 

We focus on IR emission measured by the Spitzer Space 
Telescope, emphasizing the use of 24/im to trace em- 
bedded star formation. The 24/im maps produced by 
Spitzer have very good sensitivity and resolution well 
matched to our ISM data. Many studies of nearby 
galaxies have revealed a strong empirical correlation be- 
tween 24/im emission and sites of active s tar formation 
(e,g., Ealzetti et al. 2005; Alonso-Herrero et al.l 120061 : 
[Pirez- Gonzalez et al. 2006; Prescott et al. 2007) and it 
has been used extensively to investigate the correlations 
betw e en gas and star f ormation (e.g., iKennicutt et al 
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SINGS and LVL legacy programs, Spitzer produced well- 
calibrated maps, publicly available maps for each of our 
targets. Based on these considerations, we focus on 
Spitzer data here but expect that similar studies using 
the short wavelength Herschel bands will yield comple- 
mentary results in the near future. 

2. DATA 

We focus on Esfr to complement our database of 
H I and CO maps. We therefore present measurements 
for galaxies meeting the following criteria: 1) a HERA- 
CLES CO map containing a robust detection, 2) Spitzer 
data from 3.6-160/im, and 3) inclination < 75°. This 
yields the sample of 30 disk galaxies listed in Table [TJ 
This table also gives the distance, the physical resolu- 
tion corresponding to the (limiting) 13'' angular reso- 
lution of HERACLES at that distance, inclination, po- 
sition angle, optical radius, and the source of the Ha 
data that we will use for that galaxy. We adopt redshift- 
independent distances and orientatio n from the careful 
literature co mpilation in |Kennicutt et al.l (|2011l ) wher- 
ever possible . IWalter et al.l (I2008D elsewhere, and from 
LEDA (P rugniel Heraudeaulll998f ) and NED when nei- 
ther are available. 

We will estimate Ssfr from combinations of broad- 
band infrared (IR), UV, and Ha emission with ref- 
erences to other estimates from the literature. By 
construction HERACL ES and THINGS overlap SINGS 
(jKennicutt et alJIgOOgjWA^^ et aU (j2009f ). and the 
GALEX NGS (jCil de Paz et al.l l2OO70 surveys, so Ha, 
UV, and IR are readily available for each target. 



We aim to compare tracers of the ISM and star forma- 
tion. The maps of the ISM limit the resolution of such 
comparison; the HERACLES CO maps have 13'' reso- 
lution and the naturally-weighted THINGS H I data are 
often comparable. For most of our targets, this angular 
resolution corresponds to < 1 kpc and we adopt 1 kpc as 
our working resolution. We convolve each map to have a 
symmetric gaussian beam with FWHM= Ikpc. For the 
Spitzer 24/im maps we first convert from the MIPS PSF 
to a 13'' gaussian beam using a kernel provided by K. 
Gordon (priv. comm.) then we convolve to 1 kpc. This 
effectively places our targets at a common distance but 
does not account for foreshortening along the minor axis. 
Five galaxies are too distant to convolve the ISM maps 
to 1 kpc. We mark these in Table [T] and include them 
in our analysis at the 13" angular resolution of the CO 
data. 

From these convolved maps, we generate a database of 
intensity measurements. We sample each map using an 
hexagonal grid spaced by 0.5 kpc, i.e., one half-resolution 
element. At each sampling point we measure the CO J = 
2^1, 24/im, FUV, NUV, Hi, and Ha intensities. We 
also note dust properties, estimated at coarser resolution, 
for each point. Table [2] summarizes our data set. The rest 
of this section explains how we derive the measurements 
for each point. 

2.1. HERACLES CO 

The HERA CO Lin e Extragalactic Survey (HERA- 
CLES, Lerov et al. 2009) used the Heterodyne Receiver 
Array (HERA, Schuster et al. 2004) on the IRAM 30m 
telescope to map CO J = 2 ^ 1 emission from 48 nearby 
galaxies. The HERACLES cubes cover out to r25 with 
angular resolution 13" and typical la sensitivity 20 mK 
per 5 km s~^ channel. Leroy et al. (2012, in prep.) 
present the full data set and more detail^. 

We estimate H2 mass surface density, Sh2, from CO 
J = 2 — )► 1 intensity via 



Sh2[M0Pc-2] =6.3 

(1) 

where R21 is the CO(2^1)-to-CO(1^0) line ratio and 
aco is the CO(1^0)-to-H2 conversion factor. The for- 
mula includes a factor of 1.36 to account for helium. 
We adopt a line ratio of O.TEH, and a Galactic con- 
version factor, aco = 4.4 Mq pc~^ (K km s~^) 



equivalent to Xco 



2 X lO^^cm-^ (Kkms- 



This value is intermediate amo ng recent determinations 
of the Milk y Way aco (e.g. iStrong fc Mattox "1996|; 

U WM iHever et al.l 120091 : lAbdo et al. 201(1 



Dame et al 



IHever et 
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all 120091 : ^ 
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and represents a commonly adopted "bes t single value" 
for work on ne arby disk galaxies (e.g., IWong fc Blitzl 
I2OO2I : iLerovlT al. 2008). HERACLES does span a range 

10 All HERACLES data are publicly available from 
the IRAM and NRAO web pages. Current URL 

www . cv . nrao . edu/$\sim$aleroy/HERACLES 

This line ratio is slightly lower than the 0.8 used by 
Lerov et al. ( 2009), reflecting the revised efliciency used in the data 
reduction. This value, 0.7, is the mean ratio of integrated CO J = 
2^1 HERACLES flux divided by the CO J = 1 flux mea- 
sured by lYoung et al.i |1995|), iHeliFer et al.l (l2003l ). or IKuno et al.) 
((20071 ). 
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TABLE 1 
Sample Properties 
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NGC0337 


19.3^ 


1.24 


51 


90 


1.5 


SINGS (Kennicutt et al. 2003) 


-11.43 


0.23 


NGC 0628 


7.2 


0.46 


7 


20 


4.9 


Palomar Las Campanas Atlas 


-10.84 


0.35 


NGC0925 


9.1 


0.59 


66 


287 


5.4 


Boselli & Gavazzi (2002) 


-11.10 


0.20 


NGC 2403 


3.2 


0.21 


63 


124 


7.9 


Boselli & Gavazzi (2002) 


-11.76 


0.36 


NGC 2841 


14.1 


0.91 


74 


153 


3.5 


Palomar Las Campanas Atlas 


-11.53 


0.61 


NGC 2903 


8.9 


0.57 


65 


204 


5.9 


Hoopes et al. (2001) 


-10.71 


0.56 


NGC 2976 


3.6 


0.23 


65 


335 


3.6 


LVL (Dale et al. 2009) 


d 


0.13 


NGC 3049 


19.2^ 


1.24 


58 


28 


1.0 


SINGS (Kennicutt et al. 2003) 


-11.93 


0.40 


NGC 3184 


11.8 


0.76 


16 


179 


3.7 


SINGS (Kennicutt et al. 2003) 


-11.12 


0.52 


NGC 3198 


14.1 


0.91 


72 


215 


3.2 


Palomar Las Campanas Atlas 


-11.4 


0.30 


NGC 3351 


9.3 


0.60 


41 


192 


3.6 


LVL (Dale et al. 2009) 


d 


0.62 


NGC 3521 


11.2 


0.72 


73 


340 


4.2 


SINGS (Kennicutt et al. 2003) 


-10.81 


0.57 


NGC 3627 


9.4 


0.61 


62 


173 


5.1 


SINGS (Kennicutt et al. 20031 


-10.75 


0.54 


NGC 3938 


17.9^ 


1.15 


14 


15 


1.8 


SINGS (Kennicutt et al. 2003} 


-11.25 


0.42 


NGC 4214 


2.90 


0.19 


44 


65 


3.4 


Hunter & Elmegreen ( 2004) 


-10.77 


0.16 


NGC 4254 


14.4 


0.93 


32 


55 


2.5 


GoldMine (Gavazzi et al. 2003) 


-10.89 


0.45 


NGC 4321 


14.3 


0.92 


30 


153 


3.0 


GoldMine (Gavazzi et al. 2003) 


-11.08 


0.43 


NGC 4536 


14.5 


0.94 


59 


299 


3.5 


GoldMine (Gavazzi et al. 2003) 


-11.36 


0.45 


NGC 4559 


7.0 


0.45 


65 


328 


5.2 


SINGS (Kennicutt et al. 2003) 


-10.97 


0.28 


NGC 4569 


9.86 


0.64 


66 


23 


4.6 


GoldMine (Gavazzi et al. 2003) 


-11.39 


0.99 


NGC 4579 


16.4^ 


1.06 


39 


100 


2.5 


GoldMine (Gavazzi et al. 2003) 


-11.49 


0.62 


NGC 4625 


9.3 


0.60 


47 


330 


0.7 


LVL (Dale et al. 2009) 


d 


0.55 


NGC 4725 


11.9 


0.77 


54 


36 


4.9 


Knapen et al. (2004) 


-11.42 


0.38 


NGC 4736 


4.7 


0.30 


41 


296 


3.9 


Knapen et al. (2004) 


-10.72 


0.71 


NGC 5055 


7.9 


0.51 


59 


102 


5.9 


LVL (Dale et al. 2009) 


d 


0.50 


NGC 5194 


7.9 


0.52 


20 


172 


3.9 


Boselli & Gavazzi (2002) 


-10.42 


0.60 


NGC 5457 


6.7 


4.3 


18 


39 


12.0 


Hoopes et al. (2001) 


-10.22 


0.54 


NGC 5713 


21.4^ 


1.38 


48 


11 


1.2 


SINGS (Kennicutt et al. 2003) 


-11.63 


0.55 


NGC 6946 


6.8 


0.44 


33 


243 


5.7 


SINGS (Kennicutt et al. 2003} 


-10.42 


0.54 


NGC 7331 


14.5 


0.94 


76 


168 


4.6 


SINGS (Kennicutt et al. 2003) 


-11.09 


0.61 



^ Distances adopted, in order of preference, from compilations by Kennic utt et al.l (|201in . I Walter et al.l (|2QQ8I ). and the 
NED/LEDA databases. 

^ Linear resolution corresponding to 13.3^^ angular resolution at the distance of the target. 
^ Too distant to convolve to 1 kpc resolution. Included in analysis at native (13^^) resolution. 
Flux calibration of LVL Ha maps taken to be correct. 

of metallicities (e.g.. iMoustakas et alJl2QlQ[ ) and obser- 
vational evidence suggests that the CO-to-H2 conver- 
sion factor v aries as a function of metahicity (e.g., see 
summary in iLerov etld] |2Q11[ ). However, most low- 
metahicity, high-aco systems also tend to be dominated 
by atomic gas and in this paper only the total (H 1+ H2) 
gas supply will be relevant. We therefore neglect metal- 
licity variations in aco for purposes of estimating the 
dust-to-gas ratio and Ssfr, though we consider these in 
our subsequent comparison of Sh2 and Ssfr (Leroy et 
al., in prep.). 

2.2. THINGS and Supplemental H I 

We assemble H I maps for all targets, which we use to 
mask the CO, estimate the dust-to-gas ratio, and explore 
24/im cirrus corrections. These come from THINGS 
(|Walter et al.ll2QQ8l ) and a collection of new and archival 
VLAE3 data (including our programs AL731 and AL735). 
These supplemental H I are C+D configuration maps 
with resolutions l?>"-2b" . We reduced and imaged these 
a standard way using the CASA package (see Leroy et al. 
in prep.). In a few cases the native angular resolution of 
the H I maps corresponds to a spatial resolution coarser 
than 1 kpc. In these cases we assume the H I to be 



The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 



smooth at scales smaller than the resolution. Given the 
low dynamic range in H I colu mn densities observed at 
these resolutions ( Walter et al.ir2QQ8) this approximation 
should have minimal impact on our results. We translate 
21-cm intensity into H I surface den sity assuming opti- 
cally thin emission (see references in lWalter et al.ll2QQ8f ) 
and include hehum when quoting atomic gas surface den- 
sity, Shi- Thus 

Shi [Mo pc"^] = 0.020 /m [K km s"^] (2) 

2.3. GALEX UV 

For 24 galaxies, we use NUV and FUV maps from th e 
Nearby Galaxy Survey (NGS, iGil de Paz et al.l l2007l ). 
For one galaxy, we use a map from the Medium Imaging 
Survey (MIS) and we take maps for 5 targets from the 
All-sky Imaging Survey (AIS Martin e t al.ll2005f ). We 
subtract a small background from the UV maps, de- 
termined after blanking the bright (SNR> 2), extended 
emission in the map. The magnitude of this background 
is typically ~ 5 x 10~^ MJy sr~^. This corresponds to 
^SFR ~ 4 X 10~^ M0 yr~^ kpc~^, much lower than the 
typical SsFR considered in this paper. We identify fore- 
ground stars via their U V color, by-eye inspect i on, an d 
the color-based masks of iMuhoz-Mateos et al.l (|2009bl ). 
which also blank background galaxies. We correct the 
FUV map s for th e effects of Galactic extinction follow- 
ing ^SchlegdeTaD (11998.) and .Wvder et al.. (2007. ) — see 
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2.4. SINGS & LVL IR 

We use maps of IR emission from 3.6-160/im from 
the Spitzer Infrared Nearby Galaxies Survey (SINGS, 
iKennicutt et a L"2QQ3') and the Local Volume Legacy sur- 
vey (LVL, l^le et al. 2009). We mask foreground stars 
based on UV color, by-eye inspection, and t he color- 
based masks from iMuhoz-Mateos et al] (|2009b[ ). 

2.5. Literature Ha 

We draw Ha maps from the literature for all of our 
targets. Whenever possible, we draw maps from the 
LVL survey. These have clean backgrounds and agree 
very well with previous flux measurements. For the re- 
maining targets we assembled Ha maps from a variety of 
literature sources, including the SINGS data release, the 
GOLDMine database, and surveys with the Palomar-Las 
Campanas Observatories. The photometric calibrations 
for these maps were not always available. Even when 
such calibrations were present, e.g., for the SINGS Ha 
maps, the integrated fluxes of the maps often scattered 
substantially about previous measurements (a conclusion 
verified in dependently by several of us and also implicitly 
present in IKennicutt et aI]l2009L who mix spectroscopic, 
literature, and narrow-band imaging measurements for 
estimates of SINGS galaxies). Therefore for non-LVL 
maps, we set the overall flux scale of each map by search- 
ing the literature for pinning the integrated flux of the 
Ha map to match the average literature value. Table [T] 
lists the source and adopted Ha+N II flux for each Ha. 
When several Ha maps were available for a single tar- 
get, we chose among them based on a by-eye evaluation 
of the quality of continuum subtraction and overall flat 
fielding. Because our coarse working resolution, 13'', in- 
volves heavily smoothing the maps before any analysis 
these considerations are more important than seeing or 
sensitivity. 

We process the Ha map s as follows. Whenev e r avail - 
able, we apply the masks of Muhoz-Ma teos et al.l (|2009bf ) 
to blank bright foreground stars and background galax- 
ies, which often leave residual artifacts in the continuum 
subtracted image and are in any case not associated with 
our galaxies. We also mask foreground stars identified by 
eye. During convolution to our working resolution, we fill 
in these blank regions with interpolated values. 

For non-LVL images we found it necessary to subtract 
a background determined away from the galaxy. Using 
the UV and IR 24/im images, we define a mask that en- 
compasses all bright star formation in each galaxy and 
then verify by eye that this encloses bright Ha emission. 
We subtracted the median of the map determined out- 
side this region from the whole map. In a few cases that 
displayed obvious horizontal or vertical striping, we also 
used this external region to determine and subtract a 
row- or column-wise median away from the galaxy. For 
some literature images in which the field of view was 
fairly small it was not possible to determine a robust 
median far away from all IR or UV emission. In these 
cases we subtracted the mode of a histogram of the Ha 
map after carefully blanking all bright Ha emission. Af- 
ter this processing but before any correction for Galactic 
extinction or N II contamination we renormalized the 



Ha maps so that their integrated flux matches values 
adopted from the literature. 

The Ha filter includes both Ha and [N II]. We cor- 
rect for this effect usi ng th e spectroscopic ratios from 
IKennicutt eTall (j2008D and iKenmcutt et al.J (j2009f ). If 
the [N Il]-to-Ha ratio is not available, we use the 
galaxy 's ^-band magnitud e with the scaling relation 
from .Kennicutt et al.l (|2008f ) to estimate a ratio. We cor- 
rect th^jn^£s_forth^_efe of Galactic extinction follow- 
ing [SchiiiiLiEan (|1998D . 

2.6. Dust Property Fits 

Of the Spitzer far-infrared data, only the 24/im data 
reach our working resolution but the whole IR SED pro- 
vides valuable information. We derive a variety of infor- 
mation at the coarser 40'') resolution of the 160/im 
data. Most of this inf ormation com es from fitting the 
IR SED using the Draine & Li' (20 071) models. Thi s fit- 
ting resembles that in Muhoz-Mateos et al.l ()2009a[ ) but 
we include several galaxies that they did not and derive 
dust-to-gas ratios using our new CO and H I maps. 

We carry out the dust fits using the following approach. 
First, we mask the 8/ im and 24/im i mages using the 
masks of M uhoz-Mateos et al.l ()2009b[ ). We then sub- 
tract a stellar contribution from the 8/im image using 
the 3.6/im image as a te mplate. We use a sca l ing fac - 
tor of 0.269, adopted fromlMuhoz-Mateos et all (j2009aD : 
see also iHelou et aH (j2004f ). We convolve the 8, 24, 70, 
and 160/im Spitzer images and the H I and CO maps to 
all share the PSF of Spitzer 160/im. We then construct 
radial profiles of each IR image treating our targets as 
thin disks with the orientation parameters in Table [H 
Before adopting this approach, we experimented with 
fitting SEDs to each line-of-sight but the lower signal- 
to-noise ratio and the correlation among uncertainties 
the fitted dust parameters (e.g., DGR and Umin) led to 
instability in the analysis. We found the loss of informa- 
tion from assuming azimuthal symmetry and working in 
profile to be offset by the improved S/N. 

We consider only rings with S/N> 4 at each band 
and for each IR SED we calc ulate across a gr id of 
dust emission models following iDraine et al.l (|2007[ ) and 
■Muhoz-Mateos et al.l ()2009af ) . The free parameters in the 
grid are: the intensity of the ambient radiation field Umin 
in units of the local interstellar radiation field; the frac- 
tion of dust mass illuminated by a distribution of more 
intense radiation fields (the rest is illuminated by /7min); 
^PAH, the mass fraction of dust in PAHs; and the to- 
tal dust mass. To calculate we use uncertainties esti- 
mated from the scatter in the convolved maps away from 
the galaxy or 10% of the obser ved intensity, whic hever 
is higher (this approach follows iDraine et al.ll2007D . We 
identify best-fit parameters from the minimum o^i" 
grid search. 

Based on these calculations we derive the following 
quantities for each line of sight: 1) the dust-to-gas ra- 
tio, DGR, from comparison to the convolved H I and CO 
maps; 2) the 24/im-to-total infrared (TIR) luminosity ra- 
tio — defined as uLj^/Ltir for v at A = 24/im and ^tt r 
calculated using the prescription of IDraine fc Lil ()2007[ ): 
3) the 24/im emission per unit dust mass given Umin and 
QPAU for the best-fit model. Along with t/min, these quan- 
tities will allow us to explore contamination of the 24/im 
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band and the impact of dust abundance on star forma- 
tion and map 24/im emission to total IR emission. 

We derive these IR-based quantities from radial profiles 
of maps with the ~ 40^' resolution of the 160/im data. 
For this reason, we focus on ratios like the dust-to-gas 
ratio or 24/im-to-TIR ratio that we expect to vary weakly 
within a resolution element. Under these assumptions, 
these can be extrapolated to higher resolution if part of 
the ratio is known at high resolution. 

3. UV AND Ha EMISSION 

We start with the assumption that the time-averaged 
rate of recent star formation can be inferred from the rate 
of ionizing photon production or the UV luminosity over 
part of a galaxy. Ionizing photon production is driven 
by the most massive stars, which live only a very short 
time. Assuming case B recombination, the distribution 
of ionizing photons can be traced by Ha emission af- 
ter correcting for extinction. UV emission traces mainly 
photospheric direct emission from O and B stars. 

Both approaches have a long pedigree with well estab- 
lished tradeoffs and biases. Thanks to its short age sensi- 
tivity. Ha emission directly traces the most recent gener- 
ation of star formation. However, ionizing photons may 
leak from their parent regions or be absorbed by dust, 
confusing the mapping of recombination line emission to 
local star formation. UV emission probes down to lower 
stellar masses, rendering it less sensitive to stochasticity 
and variations in the initial mass function and more sen- 
sitive to star formation at intermediate ages (a few tens 
of Myr). Horizontal branch stars and scattered light may 
contaminate UV emission, and the longer time window 
to see UV emission may bias the SFR estimate to reflect 
slightly older populations. 

3.1. Adopted Conversions of UV and Ha to SFR 

We take the relation between SFR to Ha emission from 
iCalzetti et all (|2QQ7l 

SFR [Mo yr-^] = 5.3 x IQ-"^^ Lua [erg s"^] . (3) 

A sim ilar, more recent calculation by iMurphv et al.l 
(j2Qllf ) adjusts the coefficient to 5.37 x lO'^'^. Both 
conversions derive from population synthesis modeling. 
Equation [3] adopts the default (in 2007) Starburst99 
IMF, which resembles a Kroupa (2001) IMF truncated 
at 120 M0. In surface brightness units Equation [3] is 

SsFR [Mo yr-1 kpc-'] = 634 I^a [erg s"^ sr-^] . (4) 

This calibration coeffic ient in Equation [3] is about 0.66 
times that suggested bv lKennicu ti (jl998[ ) with the differ- 
ences mainly due to the adopted IMF. Given our cu rrent 
understanding of the IMF (e.g.. iBastian et ani2010f ). the 
conversion in Equation|4]should yield more realistic Ssfr 
than the Kennicutt (1998) calibration. 

We also adopt the relat ion b etween SFR and FUV 
emission from Salim et al.l ()2007[ ). 

SFR [Mo yr-^] = 0.68 x lO-^^Lpuv [erg s"^ Hz"^] . 

(5) 

In surface brightness units this is 



SsFR [M© yr-i kpc-2] =8.1x10-2 /fuv [MJy sr"^] . 

Equation [6] comes from comparison of SED modeling 
to UV flux for lar g e s et of multiband observations. It 
adopts a Chabrier (2003) IMF. The coeffic i ent in Equa- 
tion [6] is ~ 30% lower than the Ke nnicuttI ()1998f ) value, 
even after accounting for IMF differences . It is ^ 20% 
lower than the mo re recent c a lculat ion bv lMurphy et al.l 
(j2011[ ). which like iKennicuttI (|1998[ ) considers a theoret- 
ical population th at has be en continuously star-forming 
for - 100 Myr. Salim et al.l C2007) discuss this difference. 
Given the large calibr ation data s e t and observational 
grounding, we take the iSalim et al.l (|2007[ ) FUV calibra- 
tion to be correct. From the original works, the uncer- 
tainty on the coefficients due to metallicity, IMF trunca- 
tion, and star formation history appears to be ~ 10-30%. 

3.2. Effects of Discrete Star Formation Events 

Equations [4] and [6] assume continuous star formation. 
This assumption will break down as one considers sys- 
tems with low integrated SFR, either very small galaxies 
or parts of a galaxy. For such systems, a single "star 
formation rate" becomes an inadequate approximation 
and one eventually approaches the case of a single stel- 
lar population with a discrete age. We expect a square 
kiloparsec portion of an actively star-forming galaxy to 
contain more than a single single stellar population but 
investigating this limiting case yields insight into the ef- 
fect of resolution on SFR estimates. 

Figure [T] illustrates aspects of this breakdown (see also 
Figure 6 in Genzel et al. 2010). We plot Ha and FUV 
emission (at 1500 A) as a function of time after an in- 
stantaneous burst of s tar formation. The ca lculation uses 
the Starburst99 code (jLeitherer et al.iri999l ). runs out to 
100 Myr, and adopts the default evolutionary tracks and 
IMF. 

Age Sensitivity: We use this simple simulation to quan- 
tify the age sensitivity of the two tracers, which we report 
in Table |3l We measure the falloff from peak intensity 
as a function of time, noting the time that it takes to 
reach 50% and 5% of the peak. We also derive the time 
needed to emit 50% and 95% fraction of the cumulative 
emission given off over the whole simulation. Finally we 
calculate the luminosity- weighted average time after the 
burst at which a photon is emitted, defined as 

(r) = J rLdr/ J Ldr . (7) 

Both FUV and Ha emission emit most of their light 
within a few Myr after the burst. FUV emission then has 
a long "tail" (in time) over which it continues to emit at 
a low but significant level. Meanwhile Ha drops precipi- 
tously before the burst is 10 Myr old. We re cover the ex- 
pected time sensitivity of a few Myr for Ha (jVacca et al.l 
a996';'McKee & Williams 1997) while FUV covers a wide 
range of times with a characteristic value of ~ 10-30 Myr 
and most emission gone by 65 Myr. 

Intrinsic Scatter: Luminosity varies with time after 
the burst, but the same stellar population produces the 
luminosity at all times. Considering SFR ~ M/St over 
a relatively long St ~ 100 Myr, the time- average SFR is 
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TABLE 2 

Measurements For Each Sampling Point 



Quantity 



Origin 



FUV 
NUV 



GALEX (Martin et al. 2005 : Gil de Paz et al. 2007) 
GALEX (Martin et al. 2005 : Gil de Paz et al. 20o7) 
Literature (Table [1} 

Spitzer SINGS & LVL(Kennicutt et al. 2003; Dale et al. 2009) 
IRAM 30-m HERACLES (Lerov et al. 2009, Leroy et al. in prep.) 
VLA THINGS (Walter et al. 2008) + supplemental H I 



^Dust 


Spitzer IR - 


- Draine S> 


z Li (2007) 


Dust-to-Gas Ratio 


Spitzer IR - 


- Draine S> 


z Li (2007) 


24/im-to-TIR ratio 


Spitzer IR - 


- Draine & 


z Li (2007) 




Spitzer IR - 


- Draine & 


^ Li (2007) 


24/im/EDust for f/min 


Spitzer IR - 


- Draine <^ 


z Li (2007) 
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10" 



0,01 
10* 
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10' 

Time After Burst [yr] 



10* 



TABLE 3 

Calculations for a Simple Burst Model 



Quantity 



FUV 



Fig. 1. — Behavior of a Burst of Star Formation. Starburst99 calculations for UV (blue) and Ha (red) emission from an instantaneous 
burst of star formation. The left panel shows cumulative emission as a function of time after the burst, so that one corresponds to all 
emission over the calculation. The right panel shows intensity of the population relative to intensity at 1 Myr as a function of time. Dashed 
lines in both panels indicate the 50% and 95% (left) / 5% (right) levels. 

luminosity to SFR. We calculate this scatter as an esti- 
mate of the intrinsic uncertainty in estimating an SFR 
in a regime better described by discrete events. 

We assume that a tracer is visible out to time ^visible 
and that we view it at some random time t < ^visible- 
At all t we wish to recover the same SFR. If we use a 
linear conversion of luminosity to SFR then the scatter in 
luminosity for t < ^visible is the minimum uncertainty in 
the accuracy of this Hnear conversion. We try two values 
^visible- the time at which the intensity has fallen to 5% 
of its value at 1 Myr and the time by which 95% of the 
total emission has occured. We also report the results 
of the luminosity- weighted scatter, the second moment 
in log^^o luminosity weighting by Ldr, which is a more 
natural but less intuitive quantity. 

Table [3] reports these estimates. In the limiting case 
of discrete bursts we expect a factor of ~ 2 uncertainty 
(Icr) inferring SFR from Ha and a factor of 3-4 uncer- 
tainty inferring SFR from FUV. Over a square kpc, we 



Time 

... of 50% intensity at 1 Myr 

... of 5% intensity at 1 Myr 

... of 50% cumulative emission 

... of 95% cumulative emission 

... luminosity- weighted (r) 

Intrinsic scatter^ 

... out to 5% peak emission 

... out to 95% cumulative emission 

... luminosity weighted 



3.4 Myr 
5.7 Myr 
1.7 Myr 
4.7 Myr 
2 Myr 



5.5 Myr 
27 Myr 
4.8 Myr 
65 Myr 
14 Myr 



0.36 dex 0.44 dex 
0.26 dex 0.49 dex 
0.34 dex 0.63 dex 



^ RMS scatter in log of the Ho; or FUV intensity for a 
fixed-mass burst as it evolves to the indicated point. This 
will roughly correspond to the scatter in SFR estimates as 
one isolates a single stellar population. 



the same for all times in the calculation. In this light, the 
varying luminosity implies scatter in the ability to map 
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expect to average several populations so that the scatter 
in intensity should be lower than in Table [3] by ~ V^, 
where N is the number of independent populations that 
we average (e.g., see Schruba et al. 2010). 

4. HYBRID Ha+IR AND FUV+IR SFR TRACERS 

Dust absorbs Ha and UV light, leaving only ~ 20-40% 
of the original em ission visib le in a typical spiral galaxy 
(~ 1-2 mag, e.g.. lKennicuttl [l998) and reprocessing the 
remainder into IR emission. Correcting for this extinc- 
tion represents a central challenge to SFR estimation. 

Following extensive work in the literature, we will ac- 
count for extinction by combining Ha and FUV with 
IR emission. Addin g UV and IR emission has an es- 
tablished pedigree ( tBuat|[1992|; iMeurer et aT]ll995L 119991 : 
iCortese et al.|[2QQ8[ ). while the com bination of Ha and 
IR has been explored more recently (Calzetti et al. 20071 
UCennicutt et al. 2007, 2009). Both combinations offer a 
powerful way to circumvent the problem of extinction. 

Driven by the availability and quality of Spitzer data at 
A = 24/im, we work with linear combinations of Ha and 
FUV with IR emission at this wavelength. Several sub- 
tle effects complicate the simple conversion of 24/im into 
SFR. First, the 24/im band does not capture a significant 
fraction of the total IR luminosity for ra diation fields 
< 10 ^ times the Solar Neighborhood value (jDraine fc Lil 
120071 ). Second, the exact fraction of the total IR luminos- 
ity emitted near 24/im depends on the dus t size distri- 
bution, specifically the PAH fraction in the 'Draine fc Lil 
^2007) models, and the mixture of radiation fields illu- 
minating the dust. Finally, producing 24/im does not 
require ionizing photons. Indeed, stochastic heating by 
an older stellar population may re present an imp ortant 
contribution in some regions (e.g.. Murphv et aI] [2Qlli) . 
This last point is a general concern for the use of IR 
emission to estimate SFR, one that may be somewhat 
alleviated at 24/im compared to longer wavelengths. For 
this paper the critical point is that due to these compli- 
cating factors, the use of emission to infer SFRs 
has an empirical, not theoretical, foundation. 

4.1. Definitions and Approach 

We follow recent work in the field by defining a "ref- 
erence" SFR, which we assume to be correct, and boot- 
strapping the calibration of more readily observed esti- 
mates to match this reference. We draw reference SFRs 
from literature measurements at wavelengths where dust 
is only a small concern, including; radio continuum emis- 
sion and Pag emission (|Niklas et 'anil997l : iMurphv et al.l 
120111 : iCalzetti et all I2007D . We also use UV and Ha, 
corrected for extinction based on astrophysical expec- 
tations, e.g., by the "B almer decrement" or IRX-/3 ap- 
proach (iStorev fc Hummerl Il995l: l Osterbr ^ck fc Ferlandl 
120061 : iKennicuttI ri998l : ICortese et a l. 2008). In Section 
ini we give more details on these "reference" measure- 
ments and use them to estimate the calibration of the 
24/im term in hybrid IIa+24/im or UV-h24/im tracers. 
Before doing so, we define our terms (this section), inves- 
tigate the relationships among the components of these 
hybrid tracers (Section 14. 2p . and consider the effect of 
24/im emission not associated with recent star formation 
(Section El). 

Linear combinations of Ha+24/im or UV+24/im have 
the desirable properties of working independent of scale 



and breaking apart into "obscured" and "unobscured" 
terms in an easy-to-interpret way. This behavior comes 
at the cost of accuracy, as the conversion of IR inten- 
sity to extinction may depend on User (i-e., luminosity), 
scale, or the ratio of IR-to-unobscured tracer in a com- 
plex way. For example, see Cortese et al. (2008) or the 
comparison among competing, often regime-dependent 
and non-linear, con versions of mon o chrom atic 24/im in- 
tensity to SFR in ICalzetti et all (j2010f ). We adopt 
the simple approach, focusing on linear combinations in 
which the SFR is the sum of an obscured and an unob- 
scured term. 

A priori we do not know the weight to apply to 24/im 
intensity, /24, in combination with an unobscured tracer, 
^SFR, tracer 5 to rccovcr SsFR- We solvc for this quantity, 
which we refer to as iL'tracer and define as 



^^tr 



(SsFR,ref " 5]sFR, tracer) 

0.0025 /24 



(8) 



Very simply, w is the coefficient to convert 
24/im into Ssfr in combination with a tracer 
of unobscured star formation. It has units of 
Mo yr-i kpc-2 (400 MJy sr"^)"^ 
Thus 



SSFR = S 



SFR, tracer 



^tr 



^24 



400 MJy sr"^ 



(9) 



For integrated systems, one can derive w by substituting 
total SFR for Hsfr and L24 for For ease of compar- 
ison, we cho ose the normal i zation of w so that wyia = 1 
recovers the ICalzetti et all ()2007l ) result. Equation [TOl 

It will be desirable to consider roughly how much star 
formation occurs along a line of sight without too much 
emphasis on a particular calibration. We define the 
"composite" Hsfr, (^sfr), to be the median among 
SsFR calculated in a variety of ways. The quantity de- 
pends on the adopted suite of tracers, but represents a 
useful ordinate that we will plot throughout the paper. 
Our (Ssfr) is the median of 1) from Ha + An^ = 1 mag; 
2-9J3 from Ha+24/im and FUV+24/im with and with- 
out 24/im cirrus subtraction, with only cirrus associated 
with H I subtracted, and with cirrus due to twice our 
adopted radiation field (see Section [5j). 10) FUV + 
Apuv estimated from the UV s pectral slope (the "IRX- 
p relation" Muhoz-Mat eos et a l. 2009a); 11-12) FUV + 
^Fuv estimated from the TIR/FUV rat io for you n^ and 
middle-aged populations following Cortese et al.l (2008) . 

4.1.1. Existing Work on 24lJ^ni-based ^'Hybrid" Tracers 

Ha+24um in ''H II knots'': ICalzetti et all (j2007r ) and 
iKennicutt et al.l ()2007l ) demonstrated that a linear com- 
bination of Ha and 24/im intensities recover ionizing 
fluxes inferred from Paa emission for "HII knots," bright 
regions from ~ 0.05-1.2 kpc in size, and averages over the 
central pa rts of galaxies. Assuming the Paa to trace the 
true SFR, ICalzetti et al.l (j2007D found: 

SsFR [Mo yr-1 kpc-'] = 634 I^a [erg s-^ sr-^] -flO) 

^Ua refers to the extinction (in magnitudes) of Ho; emission 
by dust. Apuv is defined analogously for FUV emission in the 
GALEX FUV band. 
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Fig. 2. — Ratios among terms in hybrid SFR tracers as a function of composite {^sfr) ■ (top left) Ratio of Ha to FUV intensity along 
a line of sight, (top right) Ratio of 24/im to Ho; emission, (bottom left) Ratio of 24/im to FUV intensity, (bottom right) Ratio of 24/im to 
TIR luminosity surface density (i^Liy) calculated at 45" resolution (~ 2.2 kpc at the median distance of our sample). The red points show 
median ratio in bins of (XIsfr) with error bars indicating la scatter. 



0.0025 /24/im [MJy sr" 



The first term is Equation HI Ttie second term attempts 
to account for Ha emission obscured by dust. Tiiis cor- 
responds to WYLa = 1 (see our Equation [9l)^ ^ 

Ha+24/J^m for Whole Galaxies: IKennic utt et al. 
([2009.) compared combinations of 24/im and Ha emission 
to Ha fluxes corrected for extinction using the Balmer 
decrement. For whole galaxies, they found = 0.68; 
for the SINGS galaxies speciflcally, they found = 
0.52. The Balmer decrements used to calibrate these re- 
sults appear uncertain , at least for the SINGS galaxies 
(jMoustakas et al.ll2010f ). We ret urn to this result b elow. 

FUV+24fim in Radial Profile: iLerov et al.l (j2008f ) pro- 
posed that FUV+24/im emission could trace the recent 
SFR for large parts of the disks of nearby galaxies (see 
also [Thilker et al...2007.) . They suggested 



SsFR [Mo yr-^ kpc-'] = 0.081 /puv [MJy sr"^] {11) 



0.0032 L 



24/im 



[MJy 



sr 



They picked the 24/im coefficient, w^vv = 1-3, to match 
other estimates in radial proflle. In proflle, this esti- 
mate matched other tracers with ^ 50% scatter extend- 
ing down to SsFR a few times 10~^ Mq yr~^ kpc~'. 

4.2. Relation Among Hybrid Tracer Components 

Figures [2] and [3] show how the components of the hybrid 
tracers relate to one another. Figure [2] shows the ratios of 
Ha-to-FUV, 24/im-to-Ha, 24/im-to-FUV, and 24/im-to- 
TIR emission as a function of (Ssfr)- Black points show 
individual kpc-resolution lines of sight. Red points plot 
the median trend with error bars indicating la scatter. 
Figure [3] shows the median fraction of the total (Ssfr) 
contributed by each term as a function of (Ssfr) (left) 
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Fig. 3. — Fractional contribution and correlation among SFR terms: (left) Median fractional contribution of different terms to XlgFR cis 
a function of composite (X^sfr): (red) Ho;, (blue) FUV, and (green) 24/im emission with w = 1 (so that the Ho; and 24/im components 
roughly sum to 1.0). (right) Rank correlation among these terms for all data above a threshold (Ssfr) as a function of that threshold. 



and the rank correlation among the terms for data above 
a given (Ssfr) as a function of that hmiting value (right). 

As (Ssfr) increases, FUV emission becomes fainter 
relative to both Ha and 24/im emission. While FUV 
has about the same magnitude as the 24/im term near 
(5^sfr) ^ 10~^ Mq yr~^ kpc~^, it contributes neg- 
ligibly, < 10%, to the overah SFR by (Ssfr) ^ 
10 ~^ M0 yr~^ kpc~^. FUV exhibits a significant cor- 
relation with Ha even up to high (Ssfr), while its cor- 
relation with 24/im diminishes quickly, becoming consis- 
tent with no correlation or a weak anti-correlation by 
(Ssfr) - 10-2 Mo yr-i kpc-^. 

The ratio of 24/im to Ha emission also increases with 
with increasing (Ssfr) though the trend is weaker than 
for the 24/im-to-FUV ratio. This is consistent with star 
formation becoming increasingly embedded in regions 
with high SsFR and agrees with trends obs erved in in- 
dividual SINGS galaxies by Prescott et al.1 (|2007[ ). Ha 
contributes ~ 40% of the total SFR over the range, 
10-^ < (Ssfr) ^ 10~^ Mo yr-^ kpc^^ equivalent to 
~ 1 magnitude of extinction. We include few lines of 
sight with (Ssfr) > 10~^ Mo yr~^ kpc~^, but in these 
regions Ha contributes a small amount to the total SFR. 

Over the whole data set, (Ssfr) > 10~^ 
Mo yr~^ kpc~^. Ha and 24/im exhibit the strongest rank 
correlation among the terms. This correlation dimin- 
ishes at high (Ssfr), with the correlation between Ha 
and FUV becoming stronger than that between Ha and 
24/im at very high (Ssfr)- However, we emphasize that 
again that we have little data at very high (Ssfr)- 

The bottom right panel of Figure [2] shows fraction 
of TIR emission emerging at 24/im increases from ^ 5 
to ~ 10% across two decades in (Ssfr)- This magni- 
tude roughly agrees with the calculations of iDraine fc Lil 
(2007) and the trend has the expected sense if dust heat- 
ing increases with (Ssfr)- The small scatter in the ra- 
tio suggests that 24/im can be used to estimate the the 
bolometric IR emission (though some of the narrowness 
in the scatter arises because we measure 24/im-to-TIR at 



the coarse resolution of the 160/im data). 

Thus Figures [2] and [3] illustrate the utility of the 24/im 
band to help trace the SFR. It correlates closely with Ha, 
traces the bolometric IR emission at higher resolution 
and sensitivity than can be achieved near the peak of 
the IR SED, and offers a way to trace the dominant, 
obscured component of Ssfr- 

The figures also show how the FUV+24/im tracer oper- 
ates. Unlike Ha+24/im, where both terms contribute sig- 
nificantly across a wide range of (Hsfr), the FUV+24/im 
works as a true hybrid, blending tracers appropriate for 
different regimes: FUV traces Ssfr at low levels where 
dust extinction is weak and 24/im traces Ssfr in dusty, 
vigorously star- forming regions. 

5. A PHYSICAL APPROACH TO 24/iM "CIRRUS" 

Weak radiation fields associated with older stellar pop- 
ulations may still excite 24/im emission, leading to a po- 
tential 24/im "cirrus." Here we define "cirrus" to mean 
infrared emission from dust heated by radiation fields not 
generated by recent star formation. Such radiation fields 
may illuminate either molecular or atomic gas, so that 
our "cirrus" can come from dust anywhere in a galaxy 
(i.e., we do not use the term to refer to exclusively high- 
latitude or atomic gas). Figures [2] and [3] suggest that 
this cirrus does not dominate the 24/im over our range 
of interest. If the cirrus did dominate we would expect 
a breakdown in the correlation between Ha and 24/im 
emission and an increase in the ratio of 24/im to other 
tracers at low (Hsfr)- Neither effect appears strong in 
the data. However, even sub-dominant cirrus may repre- 
sent an important second order correction and a physical 
estimate of contamination is important. 

Proposals to identify and subtract cirrus emission 
tend to be morphological, defining the cirrus as a 
smooth background and removing it via image process- 
ing techniq ues such as iterative median filtering (e.g., see 
Greenawal t^ejalJJlQQgMT et al.l 120001 : iDale et al.l 
i2007i : >Rahman et al.l 2010f ). Such approaches have very 
little utility at 1 kpc resolution, where the contrast be- 
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Dependence of 24fim emissivity 
on PAH mass fraction 




0.1 



1.0 

PAH Mass Fraction [%] 



i 0.1 



Dependence of 24^tm ^issivity 
on radiation field " 



0.1 1.0 10.0 

Incident Radiation Field [Mathis et al. (1983) Field] 

Fig. 4. — Dependence of 24lJinn emissivity on dust properties. 
Dust emission per unit mass at A = 24/im, 624, as a function of 
dust properties in the Draine &; Li ( 2007) models. The upper panel 
shows the weak dependence of 624 on PAH mass fraction, a measure 
of the dust size distribution. The lower panel plots the dependence 
of €24 on the radiation field illuminating the dust. Gray vertical 
bars show the value at which we fix the quantity in the other panel. 
The dotted line in the lower panel shows a linear dependence for 
comparison. 



tween star-forming regions 10-100 pc in size and the em- 
bedding medium occurs almost entirely sub-resolution. 
We therefore focus on an astrophysical cirrus estimate 
that leverages our multiwavelength data. We estimate 
the amount of dust present in each resolution element 
and the radiation field not associated with recent star 
formation. Together, these give us an estimate of poten- 
tial cirrus contamination. 

The 24/im cirrus will be the product of the amount 
of dust present, Sdust, and the emission per unit dust 
mass of dust heated by sources other than recent star 
formation, 62™®. Then 



^24 



^cirrus v /I 0\ 

— ^dust • U^j 

The dust surface density depends on the gas sur- 
face density and the dust-to-gas ratio, DGR. DGR, 
in tu rn, varies wi t h me t allicity, at least to first or- 
der (iDraine et al.l 120071 : iMuhoz-Mateos et al.l l2009al : 
iLerov et al.ll2011f ) The 24/im emissivity depends on the 
dust composition and incident radiation field. In Fig- 
ure in we p lot these dependencies as they appear in the 
iDraine fcLil (j2007[ ) dust models. The top panel shows 
the weak dependence of 634^^® on qpAU^ the PAH mass 
fraction. The lower panel shows how 624^^^ varies with 
the incident radiation fielcQ, U. To good approximation 
over our regime of interest, 0.1 < U < 10, 63™^ will 
vary linearly with the radiation field not associated with 
recent star formation. 

5.1. Dust and qpAU for Our Data 

Following lDraine Lil (|2007l) we discuss U in u nits of the local 



interstellar radiation field fISRF. [Mathis et aLIHQB S") and consider 
only changes to the intensity, not the color, of the field. 



To first order /^^^^^ depends linearly on Sdust and U 
and weakly on ^pah- We need to estimate each of these 
terms for each line of sight^ 

From our fits using the IDraine fc Lil ()2007f ) models to 
the broadband IR data, we know Sdust at the coarse 
(~ 40'') resolution of the 160/im data. To estimate Sdust 
at our finer 1 kpc resolution we assume that the DGR 
varies on large spatial scales. With this assumption, we 
can take advantage of the finer resolution of our ISM 
maps to estimate Sdust from DGR x Sgas- Not much is 
known about the behavior of the DGR on kpc scales, but 
changs: es in metallicity do tend to be weak over s uch sc ales 
(e.g., Rosolowskv et al.ll2QQ8l : iMoustakas et aIll2010D . 

Using gas to trace dust on small scales also provides 
a breakdown between dust associated with H I and dust 
associated with II2. H I is less directly associa ted with 
star formation than H2 (see references in S chruba et al.l 
I2OIID and most studie s of Milky Way cirrus focus on H I 
(e.g.. iBoulanger et al.l 11996). We may be fairly certain 
that dust associated with H I and a weak radiation field 
is not driven by star formation. We therefore experi- 
ment with cirrus subtractions that use either only dust 
associated with H I or all dust. 

Our fits using the iDraine fc Lil ()2007f ) models also yield 
^PAH, which we assume to hold sub-resolution. This term 
only weakly infiuences 624 (Figure HI top panel) so the 
approximation should have little impact. 

5.2. The Cirrus Radiation Field for Our Data 

Estimating the radiation field not associated with re- 
cent star formation represent s the most challengi ng part 
of the calculation. Following lDraine et al.l (|2007t ). we fit 
our data with a two-component radiation field model. 
One component is heated by a power law distribution 
of radiation fields extending to very high [/, presumably 
due to star formation. The other component is heated 
by a single, weaker radiation field with magnitude rep- 
resented by the free parameter Umin- In our fits, Umin 
exhibits a narrow 5-95% range U = 0.7-4.0, maximum 
^ 9, and median ^ 1.35. It displays a positive correlation 
with (Esfr) (^corr ^ 0.4) and the stellar mass surface 
density (rcorr ~ 0.4) and an anti-c orrelation with galac- 
tocentric radius (rcorr ^ 

-0.3). SeelM uhoz-Mateos et al.l 

()2009a| ) for more details. 

Umin can be interpreted as the ambient radiation field 
due to old stars but this interpretation is not unique. In 
the fit this is only the minimum radiation field illumi- 
nating the dust. A high C/min may still be driven by star 
formation so that discarding all emission associated with 
Umin may represent an overestimate of the cirrus. For in- 
stance, a square kiloparsec centered on the Sun would in- 
clude s everal OB a ssociations (Sco-Cen, Perseus, Orion; 
see Reipurth 2008a. b). These contribute substantially to 
the average ISRF so that the dominant heating for our 
local U = 1 canno t be said to com e from a truly "old" 
population (see >Mathis et al.1ll983l ). 

Because of this uncertainty, we undertake an inde- 
pendent estimate of the radiation field driving the cir- 
rus. This investigation appears in Appendix El We con- 
clude that our data suggest a cirrus driven by a ra- 
diation field t/cirrus = 0.5/7min with typical magnitude 
^cirrus ^ 0-6. Such a field would drive an equilibrium 
dust temperature ~ 1 K lower than a = 1 field, con- 
sistent with the ~ 16-19 K observed in H I and the 
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TABLE 4 
Effects of Cirrus Subtraction^ 



Description Median fdrr + 67% Range^ logio ^SFR 

Best Cirrus Estimate 0.19 (0.05-0.32) 

Only Dust with Hi 0.10 (0.03-0.17) 

Twice Cirrus 0.36 (0.10-0.58) -2.4 

^ See Section [531 

^ Median integrated cirrus fraction and 67% range by galaxy. 
^ log (^SFr) below which 95% of IR is deleted. 

outskirts of molecular clouds in the Solar Neighborhood 
(e.g., Boulanger et al. 1996; Schnee et al. 2008; Lee et al. 
I2OIII) . When we refer to "cirrus" or "cirrus subtraction" 
throughout the rest paper, we mean that we have cal- 
culated the 24/im cirrus emission associated with cirrus 
powered by t/cirms = 0.5/7min- When we refer to "double" 
or "twice" cirrus we calculated the 24/im cirrus emission 
associated with t/cirms = ^min, i-e., twice our best esti- 
mate. 

5.3. Effect of Cirrus Subtraction 

We combine our dust SED fits (to radial profiles) and 
our gas maps to estimate the 24/im cirrus everywhere 
across our sample. We assume that qpAU and the dust- 
to-gas ratio remain fixed within a radial profile and that 
^cirrus = 0.5 Umin- Combined with our higher-resolution 
gas maps this gives us a local estimate of /^^^^^ foi" each 
line of sight. When we refer to "cirrus-subtracted" or 
corrected data, this intensity has been subtracted from 
the observed 24/im intensity. 

Figure [5] displays our cirrus subtraction graphically for 
one face-on spiral. Figure [6] plots the fraction of 24/im 
emission deemed cirrus by our calculation as a function 
of composite (Ssfr)- In both plots, we illustrate the im- 
pact of our methodology by also plotting results for two 
variant cirrus subtractions: one using only dust associ- 
ated with H I and one using twice our adopted radiation 
field (see Section lS^ . Table |4] synthesizes these data into 
a few key numbers: the median and la range of cirrus 
emission subtracted by galaxy and the limiting (Ssfr) 
below which 95% of lines of sight are blanked by each 
cirrus approach. 

The cirrus subtraction suppresses faint emission at 
large radius and enhances the contrast between bright 
regions and their surroundings. Subtracting only the 
component associated with H I removes the low-lying, 
extended component but has almost no effect on the 
bright part of the galaxy where most of the star for- 
mation occurs. Setting the radiation field to twice its 
nominal value isolates only the brightest star formation, 
implying no extinction for Ha or UV emission over the 
rest of the galaxy. 

The magnitude of the 24/im cirrus identified by our 
approach resembles but somewhat exceeds the fraction 
of emission identified as powered by old stars in an anal- 
ysis o f integra ted SINGS and LVL SEDs by iLaw et al.l 
(j2011[ ). ILaw et al. (2011) found typically - 7% of 24/im 
luminosity in these galaxies, on average, to come from old 
stars. This resembles the fraction of H I-associated cirrus 
in our data but is about half the total cirrus that we find. 
Qualitatively, both studies find that dust emission from 
old stars represents a second-order, but still potentially 
important, correction in nearby disk galaxies. 



6. THE 24/iM TERM IN HYBRID SFR TRACERS 

As described in Section HI the calibration of the 24/im 
component in hybrid tracers is empirical and may vary 
with scale, target, and "unobscured" tracer. Given this 
empirical underpinning, any new approach to handling 
24/im emission requires that one re-derive, or at least 
verify, the adopted calibration w. Recall that ^^;tracer, de- 
fined in Equation [HI is the factor to linearly scale 24/im 
intensity in combination with an unobscured tracer (the 
"tracer" indicated in the subscript) to account for dust- 
obscured star formation (Section 14. ip . The lack of a 
"gold- standard" reference SFR hampers this effort, but 
we are able to draw several independent estimates of 
^SFR from the literature to check w^ia and w^w for 
our approach and sample. 



6.1. Expectation in the IR- Dominated Case 

From the typical ratio of 24/im-to-TIR emission we can 
estimate the limiting regardless of unobscured tracer, 
by considering the case where the IR term dominates 
^SFR- The median ratio of 24/im to TIR emission across 
our sample is ~ 0.07 (Figure [2]). If we equate TIR emis- 
sion with the reprocessed bol ometric light fro m an em- 
bedded, continuous starburst (jKennicuttll 1 998L corrected 
to our adopted IMF) this implies w ~ 2.4. In fact, the 
24/im-to-TIR ratio increases with increasing (Ssfr) so 
that for lines of sight where IR dominates it may be ~ 0.1 
or higher which implies a lower w ~ 1.7. Based on this 
calculation, we expect w ^ 2. It would be possible, in 
principle, to have w higher than this value for a specific 
regime but such a calibration could not successfully ex- 
tend to the case where the overwhelming majority of star 
formation is embedded and visible only through the IR. 

6.2. Ha+24iJLm 
6.2.1. Previous Work on the SINGS Sample 

Our targets heavily overlap the SINGS sample, which 
has act ed as a proving gro und for combining Ha and 
24/im. ICalzetti et al.l (j2007[ ) and iKennicutt et al.l (j2007[ ) 
determined w-^fy for Ha p e aks u sing Paa as a refer- 
ence S FR. IKennicutt et^al.l (j2007D found w^a = 1-23 
and Calzetti et al." (2007) derived wua = 1-0. This dif- 
ference presumably results from different geometry, es- 
cape fraction, and dust pro perties between M51 and the 
larp^er. Calzetti et al.l ()2007l ) sample, which included M51. 
ICalzetti et al.l ()2007[ ) found that their 24/im coefficient 
applies without modification to integrals o ver the cen- 
tral 50'' of their targets. Subsequently, Ken nicutt et al.l 
(2009) considered integrated SFRs for whole galaxies and 
found a 24/im term w^a = 0.68, and wua ^ 0.52 for the 
SINGS sample specifically. They also found significant 
s catter in wi^ry fro n i gala xy-to-galaxy. 

IKennicutt et al.l (|2009f ) interpreted their low to 
imply contamination of integrated measurements by a 
substantial cirrus component. However, they based w^a 
on Balmer decrement extinctions th at wer e later revised 
to much higher values bv .Moustakas ^eTld] ()2010l ). These 
revised estima tes should rep resent an improvement over 
those used in Kennicutt et al.i (;2009i) (J. Moustakas, priv. 
comm.). The revision affects both circumnuclear spectra 
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Fig. 5. — Illustration of the cirrus subtraction for NGC 3184- Infrared emission at A = 24/im and 1 kpc resolution. Contours run from 
0.1 to 6.4 MJy sr"-*^, stepping by factors of two. Contours > 0.8 MJy sr~^ are white. Individual panels show left to right: no subtraction, 
24/im emission after our cirrus subtraction, 24/im emission after our cirrus subtraction applied only for dust associated with H I, and 24/im 
emission after cirrus subtraction using double the recommended radiation field. 



TABLE 5 

logiQ Wlia BY Study for SINGS Galaxies 





COT 


K07 


K09 


K09+M10^ 


K09 


K09+M10^ 




Centers 


Regions 


Centers 


Centers 


Galaxies 


Galaxies 


Median \ogiQ wua^ 


0.0 


0.08 


-0.32 


0.25 


-0.32 


0.24 


Scatter (Icr) log^^Q k;hq; 


0.37 


0.15 


0.24 


0.36 


0.39 


0.52 



^ Fluxes from IKennicutt et al.l (|2009l ) and IDale et aLl 

Balmer decrements. 

^ Omitting data that yield negative or zero w. 



20071 ) with IMoustakas eFall ((20Ti 



and the radial strips used for galaxy averages EE- 

Figure [7] and Table [5] show wu g for the centra l parts 
of galaxies measured by Calzet tTet al.l (|2QQ7[ ) refer- 
encing to Paa and for both integrated measurements 
and the central parts of galaxies referenced to B almer 
decrement extinctions bv IKennicutt et al.l (j2QQ9f ). We 
also plot wi^ry from app lying the revised extinctions by 
IMoustakas et al.l (j2QlQh to the IKennicutt et al.l (2009) 
measurements. The top panel of Figure El shows a his- 
togram of all wua measurements, allowing repeats among 
galaxies using different reference SFRs. The bottom 
panel shows wua broken down by study. Table [5] re- 
ports the median a nd scatter in logio ^Ha by study. 
iCalzetti et al.l (|2007l ) did not publish their measurements 
of individual regions, but Tabl e [5] includes the M51 mea- 
surements by Kennicutt et al. (2007). However, both au- 
thors emphasize that because of their heavy image pro- 
cessing these "H II knot" calibrations should not apply 
to large parts of galaxies. 

The ensemble of measurements i n Figure [7] yie l ds me- 
dian wua = 0.9, just below the Calzetti et ajj ()2007l ) 
value. The data exhibit significant scatter, la in 
logiQi(;HQ; is ^ 0.45 dex. This scatter includes system- 
atic shifts due to choice of reference SFR and so refiects 
both an uncertainty and a true scatter. Even for a fixed 
reference SFR wua still scatters significantly (la ^ 0.15- 
0.52 dex) from galaxy to galaxy. The calibration of 



'Kennicut t et all ((2009') also work with a larger sample 
of IRAS 25/i m fluxes and Balmer decrements measured by 
[Moustakas &: Kennicu ttJ (2006). These measurements were not re- 
considered with the improved approach of Moustakas et al. ( 2010) 
but because of their higher S/N they may be less affected by 
the revision than the Moustakas et all (|2010l ') radial strip mea- 
surements (R. C. Kennicutt, priv. comm.). "Circumnuclear" 
Balme r decrements also changed signiflcantly from Kennicutt et al. 
(l2009h to Moustakas et al. (2010) and the quoted S/N for the 
IMoustakas &; Kenni cuttt i\200Gt ) data is similar to that for the 
[Moustakas et al. (2010h "Circumnuclear" spectra. 



TABLE 6 

logio 24/iM Term in the Hybrid HQ;+24/iM Tracer 



Reference SFR 



No Cirrus 



Cirrus 



Extinction-Free Estimates 

... Pa-a Centers (COT) 

... Thermal R.C. (N95,N97) 

Balmer Decrements + Ho; 
... Galaxy Centers (MIO) 
... Whole Galaxies (MIO) 
... Galaxy Centers (K09) 
... Whole Galaxies (K09) 



Other Estimates 
... Gas-Column^ 
... Ahc = 1 



- Ra 



0.05 ±0.41 
0.00 ±0.18 



0.13 ±0.34 
0.40 ±0.16 
-0.59 ±0.43 
-0.35 ±0.24 



0.61 ±0.88 
-0.06 ±0.25 



0.08 ± 0.44 
0.08 ±0.19 



0.14 ±0.27 
0.41 ±0.28 
-0.51 ±0.51 
-0.28 ±0.27 



0.58 ±0.97 
0.14 ±0.28 



Average of All Estimates 



-0.07 ±0.47 0.01 ±0.43 



^ Extinction b ased on col umn density of gas following "hy- 
brid model" in lWong Blita (^2002»V 

for large parts of SINGS galaxies thus appears very un- 
certain due to both uncertainty in the reference SFR and 
galaxy-to-galaxy scatter. 

6.2.2. wua in Our Data 

We calculate for our approach and data set. To 
do so, we use a reference SFR calculated from Ha plus an 
extinction correction. We draw the extinction correction 
from a variety of literature estimates based on , Balmer 
decrements, contrasting Fa- and Ha-emission, and con- 
trasting radio continuum and Ha emission. These are 
measured variously for either whole galaxies or central 
regions. For each extinction estimate, we calculate wua 
following EquationlHl reporting the results in Tableland 
Figure [H 

For the central parts of galaxies we use Paa-based ex- 
tinctions from .Calzetti et al.. (2007.) and Balmer decre- 
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Fig. 6. — Fraction of 24lJ^m emission identified as cirrus as a 
function of (Ssfr) from top to bottom: (green) our best-estimate 
cirrus; (red) emission only from dust associated with H I; (blue) a 
cirrus powered by a radiation field twice our best-fit value. Indi- 
vidual points show results for 1 kpc resolution elements, colored 
points show the median and la scatter after binning the data by 

(^SFR>- 



ments from 'Kennicutt et al.' (2009) and 'Moustakas et al.l 
( [2010) . Fo r full galactic di s ks wi th u se B aimer decre- 
ments from lKennicutt et al.l (|2009f ) and iMoustakas et"all 
(|2010[ ), supplement these w ith extinctions estim ated 
from the gas column following I Wong fc Blitzl ()2002L their 
"hybrid" model) and extinctions inferred from com- 
paring therm al radio continuum flux estimates from 
Niklas et all (Il995[ Il997) to our Ha maps (following 
Wong fc Blitzl I2OO2I ). For full galactic disks, we also 



benchmark against A^a = 1 mag, thou gh this is a typica l 
extinction rather than a true estimate (|Kennicuttlll998f ) . 

In the "centers" of galaxies, we work with average sur- 
face brightnesses in the central 1 kpc. For w hole galaxies 
we wo rk within r25 and we assume that the iNiklas et al.l 
(|1995[ ) fluxes come mostly from inside this area. Except 
in the case of the radio continuum, we adopt extinctions 
rather than fluxes from the literature, which somewhat 
obviates the need to precisely match apertures. Nonethe- 
less we expect our reliance on published measurements 
rather than the original data to introduce scatter into our 
results. We do not consider central measurements from 
NGC 4736 or NGC 4569 because of their Seyfert nuclei. 
We also omit Balmer decrement measures for NGC 2841 
and radio continuum estimates for NGC3351, NGC 5457, 
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Fig. 7. — Literature calibrations of the Ha + 24/im tracer in the 
SINGS sample. The factor, wuc^, to be applied to 24/im emission 
to recover the recent star formation rate in hnear combination with 
Hq; for integrals (see definition in Equation [8} over large areas in 
SINGS galaxies. From top to bottom, the plots show histograms 
of log^^o '^Hq; for (1) cill measurements, allowing repeats among 
galaxies; (2) 50 x 50^^ central regions referenced to Pao; ("COT 
Centers" Calzetti et al. 2007); (3) 20 x 20^^ central regions mea- 
sured by Kennicutt et al. (2009) and referenced to "circumnuclear" 
Balmer decrements measured by Moustakas et al. (|2010l , "K09 
Centers+MlO"); (4) the same central regions referenced to B almer 
decrement-corrected Ho: from ( "K0 9 Centers" (K ennicutt e t al.l 
12009): (5) whole galaxies measured bvlKennicutt et al. (2009) and 
corrected for extinction using the Balmer decrements measured by 
Moustakas et al. (2010, "K09+M10"). (6) whole galaxies apply- 
ing the Balmer decrements measured by Kennicutt et al.i (j200^, 
"K09"). 



and NGC 6946 because they do not return sensible ex- 
tinction estimates (they yield negative extinction). 

As in the SINGS literature, the choice of reference SFR 
and 24/im treatment affects wua, leading to ~ 0.2 dex 
systematic variation. Even for a fixed reference SFR and 
cirrus approach, we find significant scatter among galax- 
ies, typically ~ 0.3 dex (Icr). Although the 24/im offers 
a powerful, sensitive SFR tracer, its exact calibration re- 
mains uncertain. At present it appears that any adopted 
should be associated with a factor of ^ 2 uncertainty 
when applied to a specific galaxy while the average 
remains uncertain by ^ 50%. 

Figure [9] presents an alternative, more direct, visual- 
izations of our constraints on wua- We plot the Ha 
extinction drawn from the literature, Ahck, as a func- 
tion of the ratio of Ha-to-24/im intensity. These two 
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Fig. 8. — Calibration of the Ha+24fJ''m tracer: The factor, lUHa, to be applied to 24/im emission to recover the recent star formation rate 
in hnear combination with Ho; (see definition in Equation [8]) for our targets with various estimates of the true extinction. The top row 
shows the ensemble of all determinations for 24/im emission. The left column shows results for 24/im maps with no cirrus correction 
and the right panel shows i^Ha determined from 24/im maps that have been corrected for cirrus contamination. In addition to the estimates 
described in Figure [Jl we use extinction estimates based on gas column (Wong & Blitz 2002., ), thermal radio continuum ( Niklas et al..„1997., V 
and a fixed Ana = 1 mag. Galaxies do not repeat in an individual determination but repeat among estimates and in the top row. A 
vertical line indicates = 1. 



quantities directly track one another for the case of a 
fixed WYLa and the sohd curves show the expected re- 
lation for WYLa = 0.25, (bottom) 0.5,1,2, and 4 (top). 
The two panels show extinction estimates based on ra- 
dio continuum and Pa emission (top panel) and Balmer 
decrements (bottom panel) with the individual studies 
labeled. If a fixed perfectly described our data, we 
would expect to see them strung out along one of the 
solid lines in Figure [9l This view highlights the contra- 
dictions among recent studies and the large scatter in 
measurements for individual studies and galaxies. 
Treating all determinations equally, we find an aver- 
age wua ^ 0.85 with 0.47 dex scatter among individ- 
ual determinations. Taking = 1 mag and no cir- 
rus subtraction also suggests about this magnitude, 
confirming it as a reasonable starting point. This en- 



semble approach agrees with the average from the 
SINGS literature. We recommend wua ^ 0.9 as a start- 
ing point for most analyses if wua is otherwise uncon- 
strained. However, given the enormous scatter in results, 
we can not si g nifica ntly distinguish this w^a from the 
iCalzetti et al.l (120071 ) value, w^aipha = 1, and that also 
represents a reasonable assumption. 

When correcting for 24/im cirrus contamination, we 
adopt a higher = 1.3, driven by the Paa and radio 
continuum estimates. These robust, extinction- free inde- 
pendent tracers of ionizing photons suggest ^1.3, 
similar to the wu n = 1 -23 found for H II knots in M51 by 
^ennicutt et all (120071 ). Again this roughly matches the 
^Ua ^ 1 rnag case. The upper panel in Figure [9] shows 
the overall consistency of the radio continuum and Fa 
measurements with this adopted wua = 1-3, but also 
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Fig. 9. — Extinction as a Function of HoL-to-24lJ^nn Ratio: Ho; 
extinction drawn from the literature (y-axis) as a function of Ha-to- 
24/im ratio in our data (x-axis) measured in star formation surface 
density units (taking = 1 for the ratio). Sohd hues plot 

the expected relation for fixed values of iuhq: — 0-25 (bottornV 
0.5, 1.0, 2, and 4 (top); see the definition of w in Equation [8] 
The top panel shows measurements for extinction estimated from 
radio continuum and Vol measurements. The bottom panel shows 
extinctions estimated from Balmer decrements, marked by study. 
All measured HQ;-to-24/im ratios have been corrected for cirrus 
emission. The dashed gray line shows our adopted K;Ha = 1-3. The 
different studies shown here all target the same set of galaxies, so 
that figure illustrates the uncertainty in the determination of iuhck 
and the contradictory results achieved by recent studies. 

highlights the large scatter among even these "good" 
data. With better reference measurements, it should be 
possible to distinguish whether the scatter reflects the 
inadequacy of a linear hybrid to predict extinction or 
simply inconsistency among challenging measurements. 

As in the SINGS literature, the results using Balmer 
decr ement extinc tions are confusing and contradictory. 
The lKennicutt et al. (2QQ 9|) values y ield th e low-outlying 
WYicx while the Mou stakas et al.l (|201Q[ ) values yield 
among the highest w^io,. This difference holds for whole 
galaxies and galaxy centers, an so does not appear ex- 
clusively driven by weighting or S/N issues. The issues 
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Fig. 10. — Expectation for wpjjy /wiict as a function of Ha ex- 
tinction. We plot the ratio of the calibration of 24/im term in 
combination with FUV, wpuy, to the calibration of the 24/im term 
in combination with Ha, i^Ha? cis a function of the Ho; extinction 
along the line of sight. See the definition of w in Equation [8] The 
calculation assumes an extinction law with Apuv/^R = 8.24/2.33 
and a fixed stellar-to-nebular extinction ratio of Ana/^R ~ 2. 

with these measurements can be seen directly in the lower 
panel of Figure O where matched Ha-to-24/im ratios cor- 
respond to dramatically different extinction estimates. 
We incline towards interpreting this wide spread to indi- 
cate that these data do not offer a robust reference SFR . 
However, if one accepts that the iMoustakas et aT l (l2 QlQf ) 
values accurately reflect the true extinctions, the Balmer 
decrements argue for a high ^ 1.4-2.4, near the Hm- 
iting case expected for heavily embedded star formation 
(Section EJ]). 

6.3. FUV+24iJLm 

6.3.1. Expectation Relative to Ha 

Given an extinction curve and a typical stellar-to- 
nebular extinction ratio, one can r elate the FUV and 
Ha extinctions, Apuv and ^Ha- iLerov et"al] (|2QQ8l ) 
adopted an i?- band stellar-to-nebular extinction ratio o f 
A^a/AR ^ 2 (ICalzetti et al.l 119941 : iRoussel et all I2Q05D 
and an extinction law where AfuvM r = 8.24/2.33 
(ICardelH et al.l [19891 : iWvder et al] [20071 usins: their fil- 
ter definitions for R and the GALEX FUV band). They 
then solved for ^^;FUv/'^Ha as a function of extinction 
(their Equation D9), which we plot in Figure [TOl At 
high extinctions, virtually all emission is embedded and 
for FUV+24/im to match Ha+24/im then wy\jw ^ '^Ha- 
At low extinction, the ratio i^FUv/'^Ha ^ 1-8, which is 
the ratio of FUV extinction to Ha for the adopted extinc- 
tion law and stellar-to- nebular ratio. In the intermediate 
regime, near the Ahq; ^ 1 magnitude typical of our sam- 
ple, we expect t^^ruv/'^Ha ~ 1.3. 



6.3.2. 



in Our Data 



We also calculate k^fuv- For this exercise, we ex- 
pand the suite of reference SFRs and carry out the 
calculation for each line of sight with (Hsfr) > 3 x 
10~^ M© yr~^ kpc~^. Figure [11] shows histograms of 
wyw for different reference SFRs and reports the me- 
dian and scatter for each approach. The first two panels 
show wy\jy from referencing to Ha+24/im using the indi- 
cated 24/im processing and wnot] we use the same 24/im 
approach for the Ha+24/im reference and the solution 
for K^FUV- The last four panels reference to Ha with 
^Ha = 1 mag, FUV corrected using Afuv inferred from 
the UV spectral slope following .Muhoz-Mateos et al.l 
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Fig. 11. — Calibration of the FUV+24.lJim tracer: The factor 
lupuv needed to combine FUV and 24/im, calculated for individ- 
ual lines of sight (see the definition of w in Equation [8]). The top 
row shows lupuv derived referencing to HQ;+24/im with and with- 
out cirrus. The middle row shows k^fuv derived from Ha plus 
a typical 1 magnitude of extinction and FUV plus extin ction es- 
timated via the IRX-/3 relation of Munoz-Mate os et al.l ©OOa). 
The bottom row sho ws lupuv derived from FUV+TIR following 
ICortese et al.l (|2QQ8I ) for a young ("late") and intermediate age 
("mid") population, bootstrapping TIR from the 24/im. We list 
the median log^^Q lapuv and Icr scatter among lines of sight for 
each approach. A vertical line indicates wy\jy = 1- 



(|2QQ9al the so-called "IRX-/3" relation), and FUV cor- 
recte d using ^fuv inferred from the TIR/FUV ratio 
from Cortese et all (j2QQ8 ). Cortese et aL (2008 ) consider 
stellar populations with different ages. We plot results 
for their youngest population, which we call "late" (their 
r > 8 Gyr case), and an intermediate age case, which we 
can "mid" (r = 5 Gyr). 

Referencing to Ha+24/im after cirrus s ubtraction, we 
find i^ Fuv ~ 1.7. The average of the two ICortese et al.l 
()2008f ) treatments give a similar value, w^\]y 1-6. This 
i s hig her than the w^w recommended by iLerov et al.l 
(|2Q08f ) but note that we have also increased iL'Ha to 1.3 
from the w^^a = 1.0 used in Leroy et al. (2008). The 
revised wyvw is still ~ 1.3 w^ia^ in accordance with our 
expectations ( Figure flOl). Because this value can be de- 
rived from the ICortese et al.l (j2008D FUV+TIR approach 
alone, this calculation establishes the FUV+24/im as 
somewhat independent of Ha+24/im, and not a purely 
bootstrapped tracer. 

Different approaches yield w^uy from 1.3-2.1. The 
scatter in wyvv appears lower than in but this is 
mostly due to the fact that our references often already 
include 24/im. To the degree that we bootstrap the cali - 
bration from either our w^ia or the ICortese et al.l ()2008f ) 
approach, il'fuv carries the same ^ 0.2-0.3 dex uncer- 



tainty associated with those approaches. 

7. COMPARISON AMONG TRACERS AND UNCERTAINTY 

Having examined the constituent terms, contamina- 
tion, and calibration of hybrid tracers, we now gauge 
the practical impact of our choices on Hsfr- Figures [12] 
-[14] compare pairs of tracers as a function of (Ssfr)- 

Figure [12] shows the effect of our cirrus subtraction 
and compares FUV and Ha-based tracers. The bottom 
row shows that Ha and FUV based approaches agree 
well for matched treatments (see also Leroy et al. 20(i8|). 
The top row shows that the cirrus subtraction has a net 
effect of 10s of percent, but less than 100% across the 
range 10"^ < (Ssfr) < 10~^ M© yr"^ kpc^^ Below 
this range, the offsets can become much more severe as 
the Ha maps become unreliable and the cirrus correction 
has removed virtually all 24/im emission. 

Figure [13] shows the effect of pushing the cirrus 
subtraction to its limits. The most conservative ap- 
proach, removing only emission from dust associated 
with H I, has a very small impact on the net Ssfr- If 
we adopt double our nominal cirrus estimate then the 
24/im component of the tracer is largely removed below 
10~^ M0 yr~^ kpc~^, suppressing Hsfr by a factor of 
~ 2 (recall the typical ~ 1 mag). 

Figure [14] presents comparisons with "independent" 
tracers: a fixed Ahck = 1 mag ; FUV + ^fuv inferred 
from the UV spectral slope (Mu hoz-Mateos et al.ll2009aL 
"IRX-/3"); Ahq; inferred using the "hybrid" gas model of 
Wong fc Blitz (2002); and FUV combined with Apuv 
estimated from the TI R/FUV ratio for a young and 
middle-aged population (jCortese et al.l 120081 ). 

Some "tilt" can be seen in the comparisons, particu- 
larly the comparison with "independent" tracers. This 
usually has the sense that our IR-based Ssfr estimate 
recovers higher Hsfr than the comparison indicated in 
high SsFR, presumably heavily embedded regions. In 
any case magnitude of systematic differences rarely ex- 
ceeds a factor of 2 (0.3 dex) over 2 orders of magnitude 
in dynamic range. Only using gas to infer the extinction 
di verges drama tically at high (Ssfr), an issue also noted 
bv lLeroy et al.l (^2008) that presumably results from com- 
plex geometry. 

The scatter among different tracers in Figures [T2]-[l4] al- 
lows us to estimate the uncertainty in any given estimate 
of HsFR- For each line of sight in our data set, we take 
the median absolute deviation-based scatter of our en- 
semble of SsFR estimates about the mean value for that 
point. That is, we measure the scatter of Ssfr estimates 
about the mean value for each point. This folds some 
systematic uncertainty into the estimate because some 
of our approaches contradict one another. Inasmuch as 
each approach represents a reasonable method to esti- 
mate SsFR, however, this simple calculation should yield 
a good idea of how uncertain a particular Ssfr estimate 
actually is. Overah for (Ssfr) > 10~^ M© yr~^ kpc~^ 
we find a scatter of ~ 0.13 dex. We illustrate this in Fig- 
ure [151 plotting a histogram of the scatter of individual 
estimates about the median Ssfr for that line of sight. 

We have also argued that the calibration of the 24/im 
term in the hybrid tracers appears systematically uncer- 
tain by a factor of ~ 2. The 24/im contributes ~ 60% of 
the SFR on average, so this suggests a systematic uncer- 
tainty of order of 0.15-0.2 dex in the magnitude of the to- 
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Fig. 12. — Hybrid tracers with and without cirrus: Ratio (log^o) of hybrid star formation rate tracers as a function of composite (X^sfr)- 
The gray area shows where {^sfr) < 10~^ Mq yT~^ kpc~^, our lower hmit for a rehable Ssfr- Red circles show the median ratio in 
bins of (XIsfr) with error bars indicating la scatter in that bin; blue circles show the same where (XIsfr) < 10~^ Mq yT~^ kpc~^ and 
we consider SgFR unreliable. Black points show individual lines of sight. The top shows the effect of subtracting the IR cirrus on the 
integrated estimate. The bottom row compares Ha and FUV-based tracers for matched cirrus approaches 
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Fig. 13. — Effect of varying the cirrus approach: Ratio (log^^Q) of star formation rate tracers as a function of composite (XIsfr). The 
gray area shows where (XIsfr) < 10~^ Mq yr"-*^ kpc~^, our lower limit for a reliable XlgFR- Red circles show the median ratio in bins 
of (XIsfr) with error bars indicating la scatter in that bin; blue circles show the same where (XIsfr) < 10~^ Mq yr~^ kpc~-^ and we 
consider XlsFR unreliable. Black points show individual lines of sight. The left panels show the effect of subtracting only cirrus emission 
from dust associated with H I. The right panel shows the effect of setting the radiation field to double our adopted value. 
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Fig. 14. — Comparison with other tracers: Ratio (log^^o) ^^^^ formation rate tracers as a function of composite (X^sfr)- The gray area 
shows where (XIsfr) < 10~^ Mq yr~^ kpc~^, our lower hmit for a rehable Ssfr- Red circles show the median ratio in bins of (Ssfr) 
with error bars indicating la scatter in that bin; blue circles show the same where (X^sfr) < 10~^ Mq yT~^ kpc~^ and we consider Ssfr 
unreliable. Black points show individual lines of sight. Each panel compares one of our hybrid tracers to another approach. 
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Fig. 15. — Uncertainty estimated From scatter among tracers: 
Distribution of Ssfr — (^SFr), the difference between each indi- 
vidual SFR estimate in our data set and the mean XlgFR estimate 
for that hue of sight. The histogram includes all lines of sight 
with SsFR > 10~^ M0 yr"-*^ kpc~^. A typical scatter across the 
sample is 0.13 dex, which gives a rough empirical estimate of the 
uncertainty in individual Ssfr estimates. 

tal SFR. Inasmuch as we have aimed for a self-consistent 
stable of tracers, this overall calibration uncertainty ex- 
ists in addition to the point-by-point uncertainty derived 
from Figure [151 

7.1. Limiting (Ssfr) 
Our Ha maps are unreliable below (Ssfr) ^ 
10~^ M0 yr~^ kpc~^ and the details of the cirrus 
methodology dictates the 24/im distribution below this 
level. We therefore adopt Ssfr = 10~^ yr~^ kpc~^ 
as a limiting surface density for a robust estimate. The 
24/im and FUV maps recover signal well below this value, 
but their translation into Ssfr becomes more uncertain. 
The low-SsFR regime can help constrain many interest- 
ing astrophysical processes. Our general recommenda- 
tion in this regime is to carefully focus on the constituent 
parts of the hybrid tracers (i.e., FUV, Ha, and 24/im sep- 
arately) bot h as in tensities and II sfr e stimates (e.g., see 
iBigiel et in [201(1 Miruba etiD [2QlT) . FUV and Ha 
alone should provide robust lower limits to the overall 
Ssfr- 

7.2. Comparison to Integrated Values 

A key property of well-formed Hsfr estimates that 
span a large part of a galaxy is that they integrate 
to a sensible total SFR. Figure [16] shows the re- 
sults of integrating our maps, here Ha plus cirrus- 
corrected 24/im, and co mparin g them to the galaxy- 
average Ssfr of Calzetti et al.l (pQlQ). Note that the 
iCalzetti et al.l (1201 Q[) results are partially anchored to 
the ^ Kennicutt et al.l (|2009[ ) results, which we consider 
somewhat uncertain for these targets. The overall match 
appears very good given this uncertainty, with our esti- 
mates tending to be slightly higher in high SFR galaxy 
because of our higher adopted wua- 

7.3. Other Uncertainties 



Fig. 16. — Comparison to integrated SFRs: Comparison of SFRs 
derived integrating our Ssfr (a^-axis) to integrated SFRs from 
Calzetti et al. (2010) (y-axis). The solid line shows equality and 
dashed lines indicate deviations of ±50%. 

We have attempted a thorough investigation of Ssfr 
estimation on kpc scales but several important uncer- 
tainties remain. We discuss these briefly here. 

Diffuse Ionized Gas and Escape of Photons from H II 
Regions: Photons may escape from H II regions, pro- 
ducing Ha emission away from the emitting stellar pop- 
ulation. Estimates of the fraction of suc h "diffuse" Ha 
emission in galaxies are often ^ 30% (see I Rahman et aTl 
I2OIOI ). Most such emission still reflects ionizing photons 
produced by massive stars and so must be accounted in 
the SFR if kpc-sized estimates are to match integrated 
estimates. In theory this emission sho uld thus be added 
back to the likely parent region (e.g. jBlanc et al.ll2009[ ) 
and should not be entirely removed from the map. 

Given the cross section of H I to ionizing photons, the 
path length of ionizing photons near the Lyman limit will 
be extraordinarily short in the presence of any neutral 
gas. Despite this, diffusion over scales of a few hundred 
pc clearly occurs. We observe a warm ionized medium in 
our own Galaxy and diffuse ionized gas in other galaxies. 
However given that this flux must be conserved and our 
kpc resolution this systematic should represent a decid- 
edly second-order concern. 

Dust Absorption of Ionizing Photons: The comple- 
mentary concern may also confuse our measurements. 
Dust inside H II regions may absorb ionizing photons 
before they reach a hydrogen atom. Such photons will 
be missed by recombination line observations, including 
Paa, though any bolometric tracer will account for this 
light as it is reradiated by the dust. Robust estimates 
of the magnitude of this effect remain scarce, but llnoud 
(|2001l ) and ilnoue et al.i (;2001l) studied star-forming re- 
gions in the Milky Way and Local Group galaxies and 
suggested the dust may absorb as ma ny as half of the 
ionizing photons. Calculations by Do pita et al.l ()2003l ) 
also suggest plausible magnitudes of several times 10% 
with the effect strong est for ultracompact H II regions 
(see also lGroves et al.lboOSr i. 

We flnd good overall consistency between Ha-based ap- 
proaches, which will be sensitive to this effect, and FUV- 
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based approaches, which should not be affected. How- 
ever, to some degree we also calibrate the FUV-based hy- 
brids to match the Ha-based hybrids. The bottom right 
panel in Figure [M] compares Ha+24/im to FUV+TIR 
and does reveal a bias towards higher Hsfr for the bolo- 
metric FUV+TIR tracer. The bias becomes worse at 
lower (Ssfr) while we expect the opposite for dust ab- 
sorption of ionizing photons. Therefore we tend to at- 
tribute the disagreement mostly to the cirrus issue, but 
the data could certainly be consistent with non-negligible 
absorption of photons. The ideal test will be to compare 
a robust ionizing photon production map, e.g., thermal 
radio continuum emission at 1-10 mm or IR recombina- 
tion lines, to high resolution bolometric luminosity maps 
and stellar populations. 

UV Emission From Older Populations: Even a rela- 
tively evolved population may still produce some FUV 
emission, contaminating our Hsfr estimates. We es- 
timate potential contamination by measuring the ratio 
of FUV to median-filtered 3.6/im emission, an approxi- 
mate tracer of stellar mass (see appendices in lLerov et al.l 
[2008.) . in regions that appear dominated by old stel- 
lar populations. We measure ratios for the bulges 
of M 81, M 94 ("the Sombrero"), NGC 4725, and 
NGC 33510, all of which show bright stellar struc- 
tures with little evidence of current star formation. We 
also measure the ratio across our sample for all lines 
of sight with very low star formation, (Ssfr) < 5 x 
10~^ M0 yr~^ kpc~^, but significant stellar mass, > 
100 Mq pc~^. These ah suggest a ratio /puv/^* ~ 0.7- 
1.9 X 10-^ MJy sr-i (Mq pc-^)-!. 

We test the effect of correcting the FUV or Ssfr 
maps using a ratio at the upper end of this range, 
2 X 10~^ MJy sr~^ (M© pc~^)~^. Our adopted mass- 
to-light ratio is 

[Mq kpc-^] = 200 I3.6 MJy sr"^ , (13) 

similar to iLerov et al.l ()2008l), though revis ed slightly 
to match the maps of IZibetti et all (j2009f ). There- 
fore this ratio corresponds to /puv [MJy sr~^] ~ 
0.004 I3.6 [MJy sr~^] . From Equation [6] the implied con- 
tamination in units of Ssfr will be 

^SFR^^ [Mo yr-^ kpc-^] ^ 1.6 x 10"^ [M pc"^] . 

(14) 

Equation [14] implies that observations of stellar bulges 
(E* ^ 10^ M0 pc~^) that measure Ssfr only a few 
times 10""^ Mq yr~^ kpc~^ should be viewed with sus- 
picion. The contamination will be lower in the outer 
disks of galaxies (S^ ~ 10 Mq pc~^) in proportion to the 
lower stellar surface density. However, for our sample and 
working Ssfr limit the contribution of this truly old UV 
component can be safely neglected. The SFR associated 
with FUV emission from old stars appears to be only ^ 
3% of the total SFR where Ssfr > 10"^ M© yr-^ kpc^^ 
and contributes only ^ 10% of the FUV light. 

Contributions from intermediate age populations rep- 
resent a more complex problem requiring detailed stellar 
population analysis (a more complex extension of Section 
13]) that is beyond the scope of this paper. 

M 81 and M 94 are otherwise not included in this paper. 



Stochasticity and IMF Variations: At 1 kpc resolution 
we expect to average several star-forming regions in each 
element. For Ssfr > 10~^ M© yr~^ kpc~^, we expect 
> 5 X 10^ M0 formed over the last ~ 5 Myr in 
each element. This minimizes but does not completely 
eliminate concerns about sampling the stellar initial mass 
function (IMF). Changes in the IMF will directly impact 
our estimates but concrete evidence for such va riations 
remains amibiguous at best (jBastian et al.|[2QTQh . We do 
not consider either effect. 

8. DISCUSSION 

8.1. Hybrid SFR Tracers 

Hybrid tracers that combine an unobscured and an 
obscured (i.e., IR) component have many advantages. 
The two tracers complement one another, yielding an 
approach that can apply, in principle, to a wide range 
of environments. Maps of UV, Ha, and IR emission are 
widely available. IR emission directly probes reprocessed 
photons, avoiding the need for less reliable extinction 
tracers. Even given the uncertainty in the IR term cal- 
ibration, using a direct measure of reprocessed starlight 
makes for much more robust approaches than optical, 
UV, or gas-based approaches. Tracers constructed from 
linear combinations of these bands will be reasonably 
scale-independent and can be easily interpreted. Despite 
the significant uncertainties discussed in this paper, we 
argue that hybrid IR+UV or IR+Ha tracers do offer the 
best current combination of widely available data and 
robustness. 

24fJ^m Cirrus: The issue of infrared "cirrus" , emission 
due to heating by an older stellar population, remains a 
challenge for IR-based SFR estimates. Here we have sug- 
gested an approach based on physical dust models, with 
the cirrus emission mainly dependent on the dust abun- 
dance and illuminating radiation field. Several points 
emerge from our calculation. First, cirrus emission asso- 
ciated with molecular gas can make an important con- 
tribution especially in the bright parts of star-forming 
galaxies. This is emission from dust mixed with molec- 
ular gas but heated by weak radiation fields. Because 
the dust emitting in such a component is mixed with the 
clumpy molecular ISM this contribution may easily be 
missed by morphological approaches (e.g., median filter- 
ing) designed to reject a smooth component. 

Second, estimating the radiation field /7cirrus that is not 
associated with recent star formation is challenging. Do- 
ing so one must be careful to both avoid oversubtraction 
and catch all of the cirrus. We present a detailed at- 
tempt to estimate the appropriate U. Important checks 
on this quantity will come from high resolution fits to the 
dust SED in quiescent regions. Though not totally triv- 
ial to interpret, these will help inform our understanding 
of the "non-star- forming" U in galaxies. Our sample has 
been observed by the Herschel Key Program KINGFISH 
at higher resolution and we expect follow-up analyses 
to further illuminate this issue. Finally, though the IR- 
cirrus represents an important issue one must bear in 
mind that other systematics with similar magnitude per- 
sist, including the absolute calibration of the IR term in 
hybrid tracers, losses of ionizing photons due to dust or 
escape, and departures from the "continuous star forma- 
tion" approximation that underlies the whole concept of 
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star formation rates. 

Calibration and the Need For A Sample of Wide- Field, 
Extinction- Robust Ssfr Maps: We have shown that sub- 
stantial calibration uncertainties still plague the IR por- 
tion of hybrid SFR tracers. These lingering uncertain- 
ties and the other systematics we discuss highlight the 
pressing need for wide-field, extinction-robust maps of 
the distribution of recent star formation. Though ob- 
servationally expensive, such maps are indispensible to 
calibrate more widely used tracers like those we discuss 
here. Wide-field maps are a necessity because under- 
standing galaxy evolution requires Ssfr estimates that 
span whole galaxies, which in turn require understand- 
ing faint regions and characterizing diffuse emission, not 
only studying H II peaks. Perhaps the best prospect for 
a large sample of such maps is using ALMA, the EVLA, 
or the GB T to map thermal radio continuum em ission at 
1-10 mm (iNiklas et al.l llQQTHMurphv et al.|[2Qll ). Such 
maps avoid extinction concerns completely and can re- 
veal the distribution of ionizing photons in galaxies in 
exquisite detail. 

Other Uncertainties and Limiting Returns: We have 
noted several other important uncertainties: escape of 
ionizing photons from H II regions, absorption of ionizing 
photons by dust, contribution of an intermediate age stel- 
lar population, stochasticity or IMF variations, and the 
closely related age of discrete stellar populations. Several 
of these have a likely magnitude ~ 10% and remain diffi- 
cult to constrain. As an ensemble these suggest that one 
should bear an uncertainty of at least a few times 10% 
for any star formation rate estimates. Indeed at small 
scales the age-based scatter we discuss implies that the 
very concept of a star formation rate will not be precise 
beyond this level. 

8.2. Recommended Approach 

We recommend estimating Ssfr at ~kpc resolu- 
tion over the range 10~^ Mq yr~^ kpc~^ < User < 
10 ~^ M0 yr~^ kpc~^ in the following way. Use the broad- 
band dust SED along with the models of Draine & Li 
12007) to estimat e Ddnst: ^min : and a pAH- The empiri- 
cal fits of Muhoz-Mateos et al.l ()2009a[ ) may offer a help- 



ful simplification or check. Derive the 24/im emissivity 
per unit dus t mass for t/cirrus = 0.5 t/mirEH from the 
iDraine fc Lil (2007) models and subtract this emissivity 
times Sdust from the observed 24/im map (setting neg- 
ative regions to have I24 = 0). If using only Spitzer 
data, gas maps and the assumption of a weakly vary- 
ing DGR can provide a useful way to work at resolution 
higher than the Spitzer 160/im PSF. Combine the re- 
sulting 24/im intensity with an unobscured tracer, either 
Ha (Equation |4]) or FUV (Equation [6]) emission, using 
uJua = 1-3 or wfw = 1-7, respectively (Equation [9]). 

One may also attempt to correct the FUV for contami- 
nation by emission from old stars. Based on the colors of 
bulges and Iow-Ssfr, high regions we suggest correct- 
ing the FUV-related term in proportion to the old stellar 
population by Ssfr = 1.6 x 10~^ (Equation [14]) . In 
our sample we find this correction to have a small impact, 
but caution that it only accounts for old populations. 

This approach removes a physically plausible, lo- 

At very different resolution than this study, t/cirrus ~ 0-6 may 
be a safer choice. 



cally estimated cirrus component. It is quantita- 
tively and qualitatively consistent with other approaches 
over the range 10~^ Mq yr~^ kpc~^ < Ssfr ^ 
10 ~^ M0 yr~^ kpc~^. It yields maps that integrate 
into sensible integrated SFR maps. We suggest an un- 
certainty of ^ 0.15 dex in individual measurements, 
~ 0.3 dex in the calibration of the 24/im term for any 
given system, and ~ 0.15 dex for the average calibration 
of the 24/im term. 

In the Absence of Long- Wavelength IR Data: If one 
has only gas, 24/im, and Ha or FUV data one can still 
attempt a first-order cirrus estimate. Across our sample, 
24/im emissivity per unit gas mass normalized toU = 0.6 
and DGR = 0.01 is 



I24 {U ^0.6,DGR^0m) [MJy sr"^] ^ 1.2x 10"^ Sgas 

(15) 

Where the gas surface density, Sgas, combines H I and 
H2 and includes a contribution from helium. Taking a 
typical cirrus U field and DGR one can then attempt a 
somewhat less rigorous cirrus correction and proceed as 
above. 

If only H I data are present these can be used for a 
partial correction following the methodology above. In 
this case, we suggest a lower w^a = 1-0 or w^vv = 1.3. 
This will somewhat offset the lack of correction for cirrus 
associated with H2. 

Note that these suggestions assume the typical dust- 
to-gas ratio for our sample. For low-luminosity, low- 
metallicity systems the expected cirrus emission will be 
lower by a factor that is to first-order proportional to the 
metallicity. 

With No Multiwavelength Data: In the absence of 
multiwavelength data one can leverage the fact that in 
star-forming galaxies the deprojected H I surface density, 
Shi, appears relatively fiat. Across our data the median 
Shi ^ 6 Mq pc~^ with a factor of 2 scatter. In the ab- 
sence of any knowledge of the gas distribution and pro- 
vided that one is studying a "normal" z = star-forming 
galaxy, adopting Shi ^ 6 Mq pc~^ and then proceeding 
as above will allow a rough cirrus correction. 

8.3. Implication for SFR- Gas Comparisons 

Our calculations have several implications for measure- 
ments of the relationship between gas and star forma- 
tion in galaxies. The most basic implication is that we 
find broad consistency among different estimates of Ssfr 
at kpc resolution across a sample of 30 nearby galaxies. 
We demonstrate sensible scaling relations among hybrid 
tracer components, motivate a physical approach to sep- 
arate 24/im emission not associated with star formation, 
and verify the calibration of the IR component of hybrid 
tracers for our approach. 

Calibration Uncertainty and Gas Depletion Time: The 
uncertainty in the absolute calibration of the dominant 
24/im term implies a corresponding uncertainty in the ra- 
tio of gas to star formation. This quantity can be phrased 
as the gas depletion time, roep = ^gas/SFR, or its in- 
verse, the star formation efficiency SFE = SFR/M^g^^. 
In either case, our factor of ~ 1.5 uncertainty on the 
24/im term in a hybrid SFR estimate implies an uncer- 
tainty of ~ 25% on the absolute determination of the 
SFR and thus the ratio. Comparison to integrated SFRs 
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from lCalzetti et all (12010') bears out this estimate; these 
differ from our estimates only in treatment and cahbra- 
tion of the IR emission and display a median ^ 11% 
offset with ^ 30% scatter. This uncertainty will be com- 
pounded by other uncertainties in physical parameter es- 
timation: the CO-to-H2 conversion factor, H I opacity, 
stellar initial mass function, adopted star formation his- 
tory, etc. However, as this is a calibration uncertainty, 
one should still be able to carry out internally consistent 
comparisons below this level. 

Uncertainty in Scaling Relations: From intercompari- 
son of tracers we find that the systematic "tilt" among 
EsFR estimates does not exceed 0.3 dex across the range 
10~^ < (Ssfr) < 10~^ Mq yr"^ kpc"^ and is usually 
substantially lower. Here "tilt" means systematic change 
in the ratio of a pair of Ssfr estimates across this range. 
The systematic uncertainty in the index of any power law 
scaling relation involving Ssfr will be ~ tilt /Ax where 
Ax is the dynamic range of the other, non-EsFR variable 
over the range 10~^ < (Ssfr) < 10"-^ Mq yr"^ kpc~^. 
If SsFR ^ then the uncertainty in the the index n 
due to choice of Hsfr esimate will be ~ 0.15 x n. Thus 
for any n early linear relationship (e.g., Ssfr vs. I;h2, 
iBidel et al.ll2008L (201 if ), the implied uncertainty in the 
power law index due to uncertainty in Ssfr is ~ 0.15. 
For a steeper relationship lik e the Esfr, ^ found 
in the outer parts of galaxies (Bigiel e t al.ll201Q f) the im- 
plied uncertainty is ~ 0.3. The normalization of the scal- 
ing relation will have the overall calibration uncertainty 
already discussed. 

Small- Scale Scatter in Scaling Relations: High resolu- 
tion observations and simulations have begun to consider 
the scatter in the ratios of SFR tr a cers to g:as as a func- 
tion of scale CSchrub a et all 120101 : lOnodera etaP 120101 : 
J^eldmann et al. 2011). This scatter encodes key informa- 
tion on the evolution of star forming regions, the small 
scale correlation of star formation, and intrinsic varia- 
tions in the star formation efficiency. On the scale of 
individual clouds or clusters, one cannot escape the fact 
that star formation occurs in discrete events. We con- 
sidered the intrinsic scatter in Ha or FUV emission from 
the same population at different, but still young, ages. 
These suggest that as one isolates individual populations 
one will find intrinsic scatter ~ 0.3 dex in Ha intensities 
and ~ 0.45 dex in FUV intensity for otherwise identical 
regions with different ages and ~ 0.5 dex scatter. This 
scatter represents a minimum and will be compounded 
by any uncertainties in the estimation of gas mass and 
extinction or evolution of the region. The limited mea- 
surements available so far appear consistent with these 
rough calculations (see also iBlanc eTall [20091 ). 

9. SUMMARY 

We estimate the surface density of recent star forma- 
tion, SsFR, at 1 kpc resolution in the HERACLES sam- 
ple and examine the underpinnings of empirical hybrid 
tracers (Ha+24/im and FUV+24/im). We look at the 
magnitudes and correlations among constituent parts of 
the hybrid tracers and consider the effects of discrete 
events on the derivation of SFRs. Using dust models 
and mult i wavelength data, we make a physically moti- 
vated estimate of the contamination of the 24/im band by 
emission not associated with star formation. Benchmark- 
ing to reference SFRs drawn from the literature we derive 



the appropriate calibration to use 24/im with Ha or FUV 
given our cirrus approach and physical resolution. We 
then compare a wide range of Ssfr estimates and check 
the integrated SFRs for our targets. We present our rec- 
ommended approach and discuss uncertainties and im- 
plications for SFR-gas comparisons. 
We highlight the following points: 

1. Starburst99 simulations of an evolving single stellar 
population imply an intrinsic scatter of ~ 0.3 dex 
in Ha-based SFR estimates and ~ 0.5 dex in FUV 
based SFR estimates when isolating a single stellar 
population. 

2. We use the dust models of iDraine fc Lil (j2007[ ) to 
make physical estimates of the 24/im unassociated 
with star formation. These models suggest that 
emission along a line of sight depends weakly on the 
PAH mass fraction and linearly on both the dust- 
to-gas ratio, gas column density, and the radiation 
field not associated with recent star formation. 

3. By maximizing the impact of the cirrus subtrac- 
tion on low-lying emission without oversubtracting 
we find t/cirrus ^ 0.5 Umin to be an appropriate 
radiation field not associated with recent star for- 
mation. Here Umin is the fixed, pervasive radiation 
field illuminating a ll dust derived from fitting the 
IDraine fc Lil (j2007[ ) models. Our best-fit field has 
typical magnitude U = 0.6 times the solar neigh- 
borhood interstellar radiation field. The resulting 
cirrus subtraction removes ~ 20% of the 24/im, on 
average, across our sample. 

4. Subtracting only cirrus emitted by dust associated 
with H I does not offer an appreciable correction 
in the parts of galaxies where most star formation 
occurs. IR cirrus from molecular gas must repre- 
sent an important term, especially given that most 
dust is mixed with molecular gas in the inner, H2 
dominated parts of spiral galaxies. This emission 
will likely exhibit the same clumpy morphology as 
CO emission, making it difficult to estimate using 
image processing techniques. 

5. Lack of a robust, wide-field reference SFR ren- 
ders the absolute calibration and universality of the 
24/im term in hybrid tracers like Ha+24/im and 
FUV+24/im uncertain, even for the well-studied 
SINGS sample. We demonstrate this uncertainty 
in both the SINGS literature and our own measure- 
ments. 

6. Although there is substantial scatter in the avail- 
able reference SFRs, we verify the appropriate cali- 
bration for the 24/im portion of Ha- and UV-based 
hybrid tracers for our cirrus approach and 1 kpc 
resolution. Based on references to literature extinc- 
tion and SFR estimates, we recommend = 1.3 
and wfw = 1-7 for our resolution and cirrus ap- 
proach. We arrive at the FUV calibration refer- 
encing to both Ha+24/i m and the FUV+TIR ap- 
proach of iCortese et al.l (p008). The latter means 
that the FUV+24/im approach is not only derived 
from the Ha+24/im calibration. 
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7. A wide variety of UV, Ha, and IR-based ap- 
proaches to estimate Ssfr at 1 kpc resolution yield 
results with systematic disagreement of a factor of 
2 or less above Hsfr = 10~^ Mq yr~^ kpc~^. 

8. Below SsFR = 10"^ Mq yr-^ kpc-^ the Ra data 
become unreliable and our results depend strongly 
on the adopted approach to the IR cirrus. We rec- 
ommend this value as a useful lower limit for ro- 
bust SsFR estimation using hybrid tracers. Studies 
targeting lower levels should carefully consider the 
constituent terms. 

Our work also has implications for SFR-gas compar- 
isons, which have heavily leveraged these tracers. We 
emphasize the requirement that a well-constructed User 
map integrate to a sensible SFR and note: 

1. The uncertainty in the calibration of the 24/im 
term implies a corresponding uncertainty in the gas 
depletion time, Tdep = Mgas/SFR, of ^ 25%. 

2. Systematic variations among tracers are less than 
0.3 dex over the range 10~^ < (Ssfr) < 
10~^ M0 yr~^ kpc~^. Over this range, the implied 
uncertainty in the power law slope of any scaling 
relation User ^ will be < 0.15 x n. 

3. We expect an intrinsic SFR scatter of ~ 0.3 dex 
(Ha) - 0.5 dex (FUV) to emerge at small scales 
due to age effects. Observations that isolate single 
stellar populations should observe such scatter. 
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APPENDIX 

ESTIMATING THE RADIATION FIELD ASSOCIATED WITH THE CIRRUS 

Our fits to the dust properties yield tiie radiation field powering dust emission in each ring, including the best-fit 
single radiation field illuminating the bulk of the dust mass, /7min- The fits do not distinguish how much of this 
Umin is associated with recently formed stars and how much results from an older stellar population. Because of this 
uncertainty, we cannot simply identify all 24/im emission associated with t/min as cirrus. In this appendix, we attempt 
an independent estimate of the radiation field, Uciwus-, not directly associated with star formation. In the regime that 
we consider the intensity of the cirrus emission at 24/im, /24, depends linearly on the incident radiation field (Section 
[5] and Figure H]), so that /|^™ ex [/cirrus- 

Estimating [/cirrus requires that we define a "successful" cirrus subtraction, or equivalently a goodness-of-fit metric 
for a cirrus estimate. After considerable experimentation we settled on the following: a successful cirrus estimate 
removes most emission from faint regions hut does not overcorrect the map or delete emission actually associated 
with star formation. This means that we look for [/cirrus that significantly affects low-lying emission but does not 
systematically oversubtract. In practice we use five measurements of structure in the cirrus-subtracted 24/im map to 
assess the quality of our [/cirrus determination: 

1. The fraction of lines of sight for which the 24/im is set to zero, that is, completely removed, by the cirrus 
subtraction (this will include lines of sight where the 24/im has been oversubtracted) . 
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Fig. 17. — The structure of 24/im emission after cirrus subtraction as a function of (left column) adopted cirrus radiation field, t/cirrus 
and (right column) adopted scaling factor x when t/cirrus = x t/min- The top row shows the correlation of 24/im emission with Ha emission 
or Ana inferred from the ratio /24/^Ha- The second row shows the fraction of lines of sight where the cirrus removes all emission, and 
the Gini coefficient, a non-parametric measure of the unevenness of the distribution. The third row shows the derivatives of the fraction 
of pixels set to zero by the cirrus subtraction and the Gini coefficient as a function of the radiation field. The bottom row shows the 
median absolute difference between the cirrus map calculated from the indicated radiation field and the observed "low-lying" emission, 
with low-lying defined as all data less than the median. 



2. The median absolute difference between the cirrus estimate and the initial 24/im map for faint emission. We 
define "faint" emission as emission below the median I24 in the map. This measures the degree to which the 
cirrus matches the 24/im map over the faint lines of sight. 

3. The Gini coefficient (Gini lOH), a non-pa rametric measure of the unevenness of the distribution of emission (see 
discussion in Muhoz-Mateo s et al.ll2QQ9bl ). This measures the degree to which the IR emission is concentrated 
into a few bright regions after the cirrus subtraction. 

4. The rank correlation coefficient between Ha emission and cirrus-subtracted 24/im emission. In principle, this 
should measure the degree to which the morphology of the cirrus-subtracted 24/im map matches that of the Ha, 
a relatively unambiguous signature of high- mass star formation (but see complications below). 



Star Formation Rates in Nearby Galaxies 



27 



5. The rank correlation coefficient relating Ha emission and the ratio of 24/im-to-Ha emission, or equivalently 
Alia. This measures the degree to which the cirrus-subtracted map obeys the expectation that more actively 
star- forming regions exhibit more extinction (but see below). 

Calculations: For each galaxy, we carry out a cirrus subtraction for a range of test radiation fields and measure 
each of these five statistics for each test radiation field. We consider: (1) a single fixed radiation field, /7cirrus, for the 
whole sample; (2) a single fixed radiation field, t/cirrus, per galaxy; and (3) the case where the cirrus is powered by a 
scaled version of the best-fit Umin from the dust SED fits (i.e., t/cirms = ^^min; recall Umin is derived from a fit for each 
radial ring). We vary these test fields from /7cirrus = 0~5 and try scaling t/cirrus = ^^min with x from 0-4. We consider 
X > 1, which may seem counterintuitive, to allow the possibility that the cirrus is driven by a radiation field locally 
more intense than that indicated by the Umin parameter fit at coarser (40'') resolution or because of contribution from 
fields > Umin- In the end we do find x < 1 to yield the best results. 

Figure [TTl shows these measurements for all of our data, plotting each as a function of the radiation field t/cirrus 
assumed to power the cirrus (left column) or the factor x by which we scale /7min, ^cirrus = ^^min (right column). A 
single t/cirrus for the whole sample is certainly too simplistic, but the plots still give a good overview of the calculations. 
We include all data with (Ssfr) > 10~^ yr~^ kpc~^. We will find Ssfr to be unreliable at this low level but 
including these low-lying data in the cirrus calculations is important. A major goal of the cirrus subtraction will be 
to correct for faint but pervasive emission that extends down to this level. 

The left column of Figure [TTl shows that over the range Udrms ^ 0.3-1.5, the cirrus subtraction varies between 
creating no empty lines of sight at all to removing all emission associated with more than half of the lines of sight. 
Over this same range, the correlation of extinction with Ha surface brightness improves. The Gini coefficient also 
increases, implying a stronger concentration of IR emission into a few bright regions. The deviation between the 
calculated cirrus and the faint parts of the sample is minimized for Udrms ^ 0.4 (bottom left panel). The right column 
shows similar behavior scaling the fit radiation field (Umin) by ~ 0.3-1, with the match between the calculated cirrus 
and faint emission best for ~ 0.6 Umin- 

Identification of Udrrus' Figure [TTl suggests a rough magnitude of the cirrus. Low values, Udrrus ^ 0-3 will have 
little or no effect on the 24/im map and not produce a meaningful correction. High values Udrrus ^ 1 will suppress a 
very large part of the map and over subtract emission from low intensity regions. We want to pick a cirrus correction in 
the intermediate regime, near the regions of steep positive slope in the fraction of blank lines of sight or Gini coefficient 
where the cirrus subtraction has the maximum effect without oversubtraction. The derivative of the fraction of lines 
of blank lines of sight and the Gini coefficient appear in the middle row of Figure [TTl and indeed peak in the desired 
range. Even more basically, we argue that the most appropriate Uciwus will be the one that generates the minimum 
scatter between the calculated cirrus and the real map for low- level emission. The bottom row in Figure [T71 plots the 
median absolute difference between the cirrus and low-lying emission, defined as I2A below the median. 

Thus, we have several candidates for /7cirrus- 

1- ^cirrus that yiclds the maximum derivative of the fraction of blank lines of sight vs. U . 

2. /^cirrus that yields the maximum of derivative of the Gini coefficient vs. U . 

3. /^cirrus that yields the minimum absolute difference between the calculated cirrus and the observed 24/im intensity 
for faint lines of sight. 

4. U at which the fraction of blank pixels passes some fiducial threshold; 0.2 works well for our data based on the 
tests in the next section. 

5. U that yields the maximum correlation of Ha and 24/im. 

6. U that yields the maximum correlation of Ha with Ahck, traced by /24//HQ; after cirrus subtraction. 

Testing of Our [/-Fitting Metrics: We evaluate these metrics by testing their ability to recover a known radiation 
field. We generate a series of simulated 24/im data sets with cirrus powered by a range of [/cirrus and a similar set 
powered by scaled versions of Umin- To do so we: 

1. Scale the Ha data to correspond to the 24/im emission expected for Ahq; = 1 mag. This represents a plausible 
core of IR associated with star formation. 

2. Rescale the calculated cirrus for each line of sight to have a known Uciwus between and 5 and then add it to 
this "star forming" 24/im emission to produce a realistic mix of star formation and cirrus. We run a parallel set 
of calculations with t/cirrus = ^^min and x ranging from to 4. 

3. Add lognormal (multiplicative) scatter with magnitude 0.15 dex to the cirrus map to reflect uncertainty in the 
cirrus estimate. 
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Fig. 18. — Test of U -fitting metrics. Ability of our t/cirrus goodness-of-fit metrics to recover the radiation field powering the cirrus in 
model case where the true f/drrus is known. The lines show U derived from maximizing the derivative of the fraction of blanked lines 
of sight; maximizing the derivative of the Gini coefficient; minimizing the scatter between the cirrus and low-lying emission; matching 
a fixed fraction (here 20%) of blanked lines of sight; maximizing the rank correlation between Ho; and /24; and maximizing the rank 
correlation between Ha and /24/^Ha (a quantity closely related to the Ho; extinction, Ahq.)- The left panel shows the results of Monte 
Carlo calculations aimed at recovering a fixed t/cirrus^ the right panel explores the ability of our metrics to recover the scaling factor x when 
t^cirrus = xUmixi^ with C/min fit from the IR SED. A solid gray line in each panel shows equality for reference. 

We apply each metric to each simulated 24/im map and to estimate either /7cirrus or x. Figure [18] plots the ability of 
the methods to recover the input cirrus radiation field. The left panel plots results for a fixed /7cirrus5 the right panel 
for a scaled version of t/min, so that /7cirms = ^^min- 

The tests show that the derivative-based and minimum scatter metrics recover the input t/cirrus well with a small 
bias to recover higher values than the input. This small high bias reflects the fact that low lying Ha pervades much of 
the map so that forcing many lines of sight to zero actually moves some "star formation" emission into "cirrus." The 
effect is mild. Picking the Uciwns that recovers 20% completely blank lines of sight also works well, but this particular 
value of 20% had to be tuned for our data. 

The rank-correlation based metrics, Ha with either I2A or show strong biases, especially toward large 

^cirrus or X. The correlation of /24 with /hck is biased toward to recover lower t/cirrus than the input because the cirrus 
correction is noisy. As we apply progressively more cirrus correction to /24, this scatter artificially suppresses the 
correlation. The extinction correlation is biased to recover higher t/cirrus than input because the two quantities, 
and /24/^HQ;, are intrinsically anti-correlated for random /24. 

We conclude from this test that the derivative of the blank-pixel fraction, the derivative of the Gini coefficient, 
a tuned blank-pixel fraction threshold, and the minimum absolute deviation between calculated cirrus and observed 
intensity all represent useful tools to identify the radiation field driving the cirrus (though they have a mild bias to 
recover too much cirrus). We expect that the rank correlation coefficients, which fail in our data, will emerge as very 
useful for higher resolution studies but the resolution of our current data set is too coarse to take full advantage of 
them. 

Galaxy- by- Galaxy Estimates: So far we have examined the ensemble of our data, but the radiation field certainly 
varies among and within galaxies. Our metrics require a population of points to identify a best-fit t/cirrus- Therefore 
they cannot be applied point-by-point, but our individual galaxies do have enough data to derive U cirrus- 

Figure [19] shows the median /7cirrus (left) and the best-fit scaling of /7min (right) derived from our four useful goodness- 
of-fit statistics: maximizing the derivative of the blank line of sight, maximizing the derivative of the Gini coefficient, 
finding t/cirrus that minimizes the difference between 24/im and cirrus for faint lines of sight, and finding t/cirrus that 
blanks a fixed fraction of lines of sight. We plot both t/cirrus and the best-fit scaling x in Uciwus = xUmin as a function 
of the median fit radiation field Umin for that galaxy from the Draine fc Li ( 200 7) models. 

The left panel of Figure [19] reveals a median /7cirrus of 0.6 and a good correspondence between our /7cirrus and fit 
^min- The rank correlation relating the two is ~ 0.9, but the derived t/cirrus are smaller than the fit Umin by a factor 
of two on average (dashed line). The exceptions are the two bright systems NGC 337 and NGC 7331 where both the 
dust fitting and our metrics suggest a very high diffuse field U ~ 4-5. These targets are interesting points of follow-up 
but appear to be special cases, as we expect the very high field to be partially due to resolution and geometric effects. 
We do not follow up further here, but note these galaxies as excellent targets for Herschel study. 

The right panel of Figure [19] shows the best-fit scaling x for the case where t/cirms = ^^min, again as a function of 
the fit field Umin- The median scaling is x ~ 0.6 with about a factor of two scatter. There is certainly some correlation 
between host galaxy and best scaling and again NGC 337 and NGC 7331 appear as outliers, though here with very 
large error bars, meaning large scatter in x determined by different methods. This median scaling in the right panel 
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Fig. 19. — (left) Calculated t/cirrus vs. f/min from fitting dust models to the IR SED. We show the median f/cirrus (y-axis) estimated for 
individual galaxies (points) using our four reliable methods. We plot Udrrus as a function of each galaxy's average f/min, the fixed radiation 
yielded from fitting the Draine & Li (2007) models to the IR SED. t/cirrus tracks f/min well, but tends to be lower than f/min by a factor 
of ~ 2. The dashed line shows t/cirrus = 0.5 t/min, which is a good description of the data. The dash-dotted line show t/drrus ^ 0-6, the 
median value. The right panel shows an analogous determination of the best-fit factor x considering ?7cirrus = x t/min- The dashed line 
shows the median x ~ 0.6. NGC 337 and NGC 7331 appear as outliers in both panels. 

and the typical ratio of /7cirrus to to Umin in the left panel agree with one another fairly well, both suggesting a cirrus 
radiation field about half of the fit Umin- 

The good correspondence between our calculated Udrrus and the fit t/min suggests that we can use a scaled version 
of t/min for local estimates of the cirrus. This is desirable, as it allows us to capture local variations in the radiation 
field driving the cirrus. The correspondence of t/cirrus and Umin also reinforces that our calculations have returned a 
meaningful result, inasmuch as they return information similar to an SED fit. 

Conclusions From Radiation Field Calculations: We define and test a series of goodness-of-fit metrics that 
attempt to identify the radiation field powering the cirrus at 24/im. We require that subtracting a cirrus driven by 
this radiation field significantly affect the faint emission in the map but not dramatically oversubtract emission. 

For our data set, the cirrus correction should employ Udrrus ^ 0.5/7min, a median Udrrus ^ 0.6-0.7. There is a 
factor of ~ 2 scatter in these estimates from galaxy to galaxy with two notable outliers. In practice this lowers the 
median U used in our calculation from ~ 1.4 times the Solar Neighborhood ISRF on average to ~ 0.6. This appears 
reasonable; as noted, the Solar Neighbornood ISRF likely includes significant contributions from nearby young stars. 
When we refer to "cirrus" throughout the rest of the paper we will adopt Udrrus = 0.5 x Umin- To test the impact 
of this conclusion, we will also occasionally refer to "twice" or "double" cirrus, which indicates that we have taken 
^cirrus = 1-0 X Umin, l-G., doublc our bcst estimate of the appropriate Udrrus value. 

The critical points are that the rough magnitude of the cirrus is slightly below the Solar Neighborhood ISRF and 
does appear to track the radiation field estimated from dust fits. The specific scaling Udrrus ^ 0.5 Umin almost certainly 
results from the interaction of resolutions for which we carry out our dust fits (40'') and cirrus calculations (1 kpc~ 13'') 
and we do not expect it to be a general result . For a recent demonstration of th e shortcomings of si mple power law 
distributions of U to yield scale-independent dust-fitting results see the work by iGalliano et al.l ()201L) on the Large 
Magellanic Cloud. The magnitude, Udrrus ^ 0.6, on the other hand, seems plausible and may represent a reasonable 
conservative starting point (as might U ~ 1). In any case, our understandin g of IR cirrus emissio n will benefit greatly 
from achieving good IR SED coverage at very high spatial resolution (e.g.. iLawton et aDl2010f ). This will allow the 
kind of dust-SED based analysis that we apply here to be undertaken on maps with the ability to cleanly distinguish 
quiescent and star- forming parts of a galaxy. In this regime, we expect the rank correlation statistics that we test and 
reject {I24 with Ha morphology) and simple estimates of the total U in quiescent regions to become powerful tools. 

EsFR ESTIMATES FOR HERACLES DISK GALAXIES 

In this appendix we plot Ssfr estimates based on our data. Figures [2Ql - [24l show each line of sight considered in 
this paper as a point color-coded by Ssfr- As described in Section [21 these data are sampled on a hexagonal grid 
with spacing 0.5 kpc except for a few distant cases noted Table [U which are sampled with 6.5" (half-beam) spacing. 
Gray points indicate where the Ha+24/im based Ssfr < 10~^ Mq yr~^ kpc~^, our limit for a robust measurement 
(Section [7]). From left to right panels show User estimated from IIa+24/im, FUV+24/im, FUV only. Ha only, and 
24/im only. In all panels we have applied our recommended cirrus correction (Section [5]) to the 24/im data and adopt 
our recommended wua = 1.3 and ^^;FUV = 1-7. The 24/im-only panel uses w = 1.3. 
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Fig. 20. — 5]sFR estimates for the individual lines of sight studied throughout this paper. The axes show displacement from the galaxy 
center in kpc. We color code points by log^^Q XIsfr estimated following the recommendations in the main text. We fix the scale for each 
galaxy to allow comparison of different components but the color scale varies from target to target according to the dynamic range in the 
maps. Points where Ssfr estimated from HQ;+24/im < 10~^ Mq yT~^ kpc~^ are color coded gray in all panels. 
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Fig. 22. As Figure f20l 
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